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SAUDAÇÃO 


J. N. FERREIRA DIAS J.º% 


Professor do 1.8. 7T, 


Passou meio século sobre a fundação do Instituto Superior Técnico. Não pertenci aos primeiros 
cursos que o frequentaram e já não ouvi as aulas dadas em francês por alguns professores estran- 
geiros, enquanto não aprendiam a nossa língua em termos de se fazerem entender. 

Matriculei-me em 1918, mas ainda senti o quadro das primeiras horas. Andei pelos bancos 
dos velhos edifícios do Conde Barão ; ainda tive como director o Prof. Bensaúde, inspirador da criação 
do Instituto; ainda tive Tomás Bordalo Pinheiro como professor de Desenho; ainda obedeci à figura 
patriarcal do Tavares, chefe do pessoal menor. E ao dizer que obedeci, faço-o sem afectação, com 
a mais perfeita humildade; hoje, certa mocidade da vanguarda acha impróprio obedecer a alguém. 

Vi ainda, sobre o pano de fundo dos barracões de chapa ondulada, coisas que já não se podem 
tornar a ver: a gravata branca do Prof. Mira Fernandes, a barba mefistofélica do Prof. Lepierre, 
a farda do general ferrugento Gonçalves. Ouvi ainda, nas velhas aulas desconfortáveis, sons que 
já não se podem ouvir: o português afrancesado do Prof. Fesch, as palavras em surdina, macias 
como um zéfiro, do Prof. Santos Viegas, a voz de clarim do severo mas gentil Prof. Borges Sequeira, 
que dava notas negativas e concedia aos alunos, depois da leitura do ponto de exame, «um quarto 
de hora para adoecerem». 

E recordo algumas das tropelias que fiz: misturava batatas com os modelos de cristalografia 
(os «batatoides») nas aulas práticas de Mineralogia; dava saltos de barreira sobre a mesa do mestre, 
na aula onde ensinavam os Prof. Droz e Adrião Sequeira, aos quais peço desculpa do desacato; 
de parceria com o Mariz Simões, algumas vezes resolvi pontos de exame de Electricidade a troco 
de paus de chocolate (1). 

Mas este comércio e estes cometimentos desportivos, talvez irregulares, faziamo-los, eu e os 
do meu tempo, com a irreverência que a idade pede e desculpa; mas fazíamo-los com o coração 
puro, sem aceitarmos tráfego com ideias de desrespeito às autoridades, à escola ou à Nação. 
Eram princípios que não discutíamos; eramos apenas rapazes que brincavam — e que também estu- 


davam. 


O Instituto Superior Técnico nasceu numa hora oportuna, porque tem sido, através dos seus 
diplomados, um elemento de primeira linha no desenvolvimento económico nacional. 


(1) Depois de escritas estas palavras, morreu o Prof. Droz; curvo-me perante a sua memória. 

* Destinava-se esta «Saudação» à abertura do novo número especial comemorativo do cinquentenário do 
Instituto Superior Técnico. Devido às perturbações então havidas na vida associativa, tal edição não foi possível, 
nem seria, de momento, oportuna. 
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Quando, à volta de 1930, a frequência de alguns cursos caiu a números digitos, porque 
as actividades portuguesas não absorviam todos os diplomados, muitos responsáveis acharam que 
a escola tinha especializações a mais para o nível do País, e poucos reflectiram que era o País que 
tinha desenvolvimento a menos para as necessidades da época. 

Afinal, venceram estes; está á vista a prova de que seguiam no bom caminho. Mais do 
que isso, verifica-se que, se a escola não cumpre, é apenas porque fornece ao País diplomados 
a menos. 

A carência de técnicos para satisfazer as necessidades crescentes da Metrópole e Ultramar 
cria já hoje numerosas dificuldades; dentro de algum tempo criará impossibilidades. Este é o pro- 
blema mais grave que a premência do crescimento nacional põe à reflexão das escolas de todos os 
graus e de todos os ramos; é preciso que elas encontrem maneira de intensificar a formação das 
novas gerações e de permitir mais largamente o acesso ao ensino superior de todos aqueles que nos 
graus mais baixos se revelem prometedores. 

Com ressalva da sua parcela de responsabilidade nesta rarefacção de profissionais, o Instituto 
Superior Técnico merece — e com ele todas as escolas — o reconhecimento nacional pelo saber que 
difunde. O empirismo e o autodidactismo estão em decadência ; cada vez é mais importante o papel 
do ensino na formação dos dirigentes; cada vez são precisos mais dirigentes em cada actividade ; 
cada vez é preciso criar mais actividades. 

É-me grato reconhecer à minha escola a posição de órgão-chave nesta progressão rapidamente 
crescente, 
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C. D. U. 517.513:512.831 


Formas Matriciais em Problemas de Extremos 


por A, GOUVÊA PORTELA 


PROF. DO |. S. T. 


o — Introdução 


1 


Seja u=u (x!,...x"), uma função contínua derivável até à ordem s qualquer, em relação 
às variáveis x! .... x” numa dada região do espaço em torno de um Ponto P dado. As variáveis 
x!.... x”? são medidas num referencial ortogonal. 

Admitamos que o desenvolvimento em série de Taylor é possível e que se pode escrever : 


A PO, + dn J)=u (ocrnho )= Au Stu + O Mu. 


Sendo: 
1) iron DE] 5 [0088 com e É 
(nm) (n,1) E ad 
A u= [0x]. [uj]. [0x] o Ra crnda 
a) (1,n) (nn) (n,1) O | uj= ms 
dx! dx! 
Au = [[3x] Lug]. [ôxi]]. [9x8] TA. “a 
3) (1, n) (n, n) (n, 1) com | o? u 
ni Ata Ui me 
(1, n) (n,1 di dx! dx) d xk 
Su = [9x] [[9 x] o Liguei) o [3x']| . [9xk] TS q, 
| 4) (1, n) (n, n) (n, 1) (n, 1) com | d+ u 
—— DIA e DE DX 
(1,n) n,n) dxt dx dxé dx 


Tem interesse em substituir as formas 3), 4), etc., por outras que sejam mais fáceis ao trata- 
mento matricial e a apresentação destas novas formas constitui o objectivo do presente artigo 


I — Formas 4º e 4“ 


As formas At e 4? estão perfeitamente adaptadas ao tratamento matricial, por isso começa- 
remos por tratar as formas 4º e 4º e no capítulo seguinte efectuaremos a sua generalização. 


1,1) 4ºU pode escrever-se ainda sob a seguinte forma: 


NU = [05]. [Ux ] LEE GE] 


(L,n) (n,n?) (n2,1) 
Sendo | Us |=| [Usa] [U;3]  [Uy.n] 
| | | (n,n?) 
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| Lu “g 


e [0x dx] =| 0x. dx 


n 
dx .9xn 
Ox dx kn 


Isto é, [U;;x ] é uma matriz rectangular com n linhas e nº colunas. 
e [xi dx*] é um vector coluna com nº linhas com elementos que resultam do produto de dois 
deslocamentos (9x). 0xk). 


1,2) 4*U pode escrever-se sob a forma seguinte: 


MU = [050.090]. [Us . ml. [0% « 0x*] 
(1,n?) (n?, nº?) (n, 1) 


 EOBRA s sm [Ugan] 


Sendo: 


[U; . u] = 
e ADygnllo.<:LUgan].! 


(n?,n?) 


isto é, [U;;, w |] é uma matriz quadrada de nº linhas e nº? colunas. 
Tanto a forma 1,1) como 1,2) são fáceis de verificar por simples efectivação das operações indicadas. 


II) Generalização 
É imediata a generalização para derivadas da ordem par 2r ou ímpar 2r41 sendo r=0,... n 


2,1) Formas pares A?" u 


H 


| 2" | 
a É cas e, E a 
seua) côx =" | | 7x 


1 1 - 
(1,n") (n",n") (n",1) 
Sendo: 
2r 
pe du 
DRE caca —— E de= Ta 
q dx .... dx9 
—— -— a 
Z2T 
e 
3 x! xt.) 
E n 
O e E O 
tóxi| = | 
1 | 0 x! xt. Ixn 
A 9 2 
(1 nº) xe. 
E q om RÉ Ê n 
0 x! O 0x! | nº 
9 x 3 x! 5 x1 
n 
dx"... dx! dx 
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A matriz [U;... q] toma a forma seguinte: 


(u”,nf) 
| = n n n 
e O e O 
2 ir € LÁ or qo E E LDA, ars ME vs + “Un. in.u“Umn.an.an 
* * . o . EE . 

5 PS tu Uu .... mm . . d Hd * . - 

. : . n n “ n º : 
Uu..mi Uu..nn, Um tm E UM o BE cos BM, “Uin.ain.ni“Uímn.n.nn 

n' ; : 
n “ . 

Uai..Un Um..mn | “Ume..nn.a “ Um.nn..m “Unn...nn.. 4 Unn..nn.n 

n nje: n : 
Un:..nmi Uni... nn - -Uniy..nn..ni Un»..nn..nn -Unn..nn..m Uun..nn..nn 

mm e 
n n 
— -— — an? 
n' 


2,2) Formas ímpares A" +! uq 


r t d+ 2+7+1 
setiu =| cõx'| Er] É poi | 
54 : | 


(1, na") (nº, nº+!) (HH 7) 
Qr+1 2r+1 
PRE a d u 
Sendo Ui...q = 2 


dx! Md MM Me. dx9 


A A O 


2r-+1 


A matriz rectangular [uj....qgl toma a forma explícita seguinte : 


(nº, nº!) 
nrH 
—— — = ——— — 
pe er+ 2r+41 ' ? - | 
-——— — —— —— ——— — — . 

Us, ... t paid U,,in... In ] | " ” U,,4l ... 1 U,,in In 

| NE 7 

Ui, ... dn Ui,in...mm . Un, ... nl U,n,in...m 

——— nn é iá — DD ————— um e a —* 

nº E e Nr 

: nf hd * “ “ “ . . o n" 

Uj nt... 41 U,,nn... In : Unt.co Ml Un nn... dn 

: n 4 

U,, ou cc nf Ui, “+ ] : u, nm «co nl U, no «nn 

— o ee— a ——. 
e nº R i J 
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3) Aplicação das formas referidas em 2) 


Tem interesse aplicar as formas atrás apresentadas a casos particulares. 
O mais simples corresponde a n = 2. 


9 
== 0) O u, dx! + u de= 2 w0x 


o x? 
92 


t=1) LSu=[ix, quina. Des de | [5] = E uwjixiôx 
gi ij=1 
Us U; 9X" 
Su = [5x!,0x)]. fm Ujz Us ni | «| Oxt dxl 
Um Uizz Us Uso» 3 xl 0d x? = E uixôox dx) 0 xk 
0x? 0 x! 
ria | 
r=2) Au=[ixtix!t,ox!ix?,0x)º 0x), 0x2 3x)]. | um Um . Us Us 0 xt dx! 
Ui Uia O May Migas | 2S 0 xt 0 xº | = 
Up Usa * o Us Uzao 0 x? 0 x! 
| Us Uzços Uz22 Usos | [9 es 


= ujki.Ixidxl dxk dx 
iJkI 


13 

4) Condição suficiente para que uma forma par seja definida positiva (ou negativa) ou semidefinida positiva 
(ou negativa) 

Seja dada uma forma par: 


r tt o r di 
seru=| 23x |. |uca || 7ês | 
| 


(1,n”) (n",n”) (n"',1) 


Como a matriz [uj...q] é simétrica por construção, é sempre possível encontrar uma 
r r 
(n",n”) 


matriz ortogonal T(n",n") tal que: 


Ê T. Ti=E Sendo E uma matriz unidade. 


T. lu...g]). Tt = [D] = [Di] 
nm) (ué,n) (nf,m') Quo) dn" ,m?) 


Sendo D uma matriz diagonal. 


Então será: 


ci t 
seu=| 23x |. TT |Ua). TUT.| sós] 
— — — + —— ] 
=— 1 rã 


== | 
“ r h 
Se designarmos por =! ao elemento da linha i do vector | + Ô d | 
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e por a! um elemento da matriz 1 


aeu=| Esto | T| Ui «a. TH| > stat | 
=[3 def |. JANE El 


(1,n") 


2 
= 2 Di; . (2 tal) 
j i 
r i 2 
= & Dj. (=2x") 1) 
J 1 


e para que A” U seja positivo (ou negativo) qualquer que sejam 0x! terá de ser Dj > 0 


a - 
(j=1,...n"), porque o segundo factor t ( 7 0 “) 1) é sempre positivo (para nem todos os 
2 


Será ainda 


0 x nulos) por estar elevado ao quadrado. 


Se Dj; > O (ou < 0) a positividade (ou negatividade) será semidefinido. 
Se Dj > 0 (ou << 0) a positividade (ou negatividade) será definida. 


Está assim generalizado às formas pares superiores à quadradada, as condições suficientes de 
positividade ou (negatividade). 

O interesse da forma é evidente, uma vez que a diagonalização de uma matriz mesmo para 
ordem de 100 — 200 é empresa fácil com os calculadores modernos. 

Assim para u=4 e r=3 aordemde 4º"u é n'= 4 =464, 

Dificilmente a estacionaridade e estabilidade de uma função da física ou da engenharia se 


decide na derivada 2X r= 6. 


Nota: porque a matriz [uj...q] é simétrica, a hipótese de os Dj; poderem ser positivos 
e negativos é impossível. 

Com efeito é fácil de demonstrar que as matrizes simétricas são semidenidasfi ou definidas 
mas nunca indefinidas. 

Se a matriz [uj...q) não estiver numa forma simétrica, poderá proceder-se à sua diagona- 
lização e se os Dj; > 0 (ou < 0) então trata-se duma forma simétrica. 

Se tal não se verificar, ou a função u não é contínua ou o referencial não é ortogonal ou 
simplesmente se cometeu um erro na operação da diagonalização. 


SUMÁRIO 


Formas matriciais de ordem superior à 2.2 e sua utilização em problemas de extremos. 
Estabelecimento de condições de suficiência. 


SUMMARY 
Matricial forms with orders above the second and utilisation in problems of extremes. 


Establissement of suficient conditions. 
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NOTAS INFORMATIVAS G. D. U. 624.314.5/94 


Elementos sobre a produção e o consumo de energia 
na rede eléctrica nacional 


— Elementos extraídos das estatísticas mensais do Repartidor Nacional de € argas (RN. C.) 


Nora: As produções e os consumos das empresas do R, N. O. representam 
cerca de 94 º/, dos totais do Pais, 


I— Breve nota mensal JANEIRO 


No aspecto hidrológico, o mês de Janeiro apre- 
sentou-se, no conjunto, muito acima da média. Efec- 
tuou-se uma exportação para Espanha de excedentes 
de energia da rede portuguesa. Os diagramas semanais 
de fornecimento são acordados entre o R.N.C. e a 
UNESA. 


II — Elementos gerais (GWh) 


) 
N al 
a = 
a) Mensais - nas 
E Su 
Variação 
| 1962 | 1969 | "ie | EE Ea 
bh 4 E + 
EEE ARE CE ag = B 
Produção hidráulica (Ph) .. | 9228) 8159|-+ 17 | [5 EE 
Produção térmica (Pr)... ...| 00/ 58] — dede 
Produção total (PT)... ...| 322,8/ 3817|+ 18 
Exportações (Ex)... ...... 0,0) 16,8] — Dom | 
Importações (1) +... ...0.0. 0,0 0,U O 
Produção para con- | ' 
sumos não perman. (Peop) (')| 608 od t|— 1 
Produção para con- RE: IV—E i 
sumos permanentes (Pep) (2) | 2620 | 305,2 | +16,5(3) nergia armazenada nas principais albufeiras 
Coeficiente de hidraulicidade | 19L| 1,59 em 
rr NE SEE | | No fim do mês 
na . Gis : Albufeiras ; ion nie ndão 
| i ” a . indi 4.7. : 
(3) À profinção: paia commons: prembanias (bar) 6 eirmbiado GWh | % (3) 
pela seguinte expressão; Pc = Pr— Penp ali 1 — Ex EAR E SE 
3) O t Il d duçã r rmanen- 
Ko, Cuide qm (cóiA é URTGRDCIA dos deh Ud a ico apelo, é Paradela. + . ce. 99,6 44,8 
respectivamente de 13,7 lo» Vida NOVE s breno a e 98,1 16,6 
Salamonde , . cc. cc... 26,3 09,3 
II — Diagramas de carga dos dias característicos CAMGRAS mma pm ova 31,6 5d, 
a Cabril: . cocos e ww) BOT 90,5 
” am feira : Castelo do Bode o. ss o wtá 149,6 91,8 
17-1-62 | 16-1-63 Guilhofrei . +... E 41 49,4 
a oe DR Lagoa Comprida + . +... 24,3 (2) 13,8 
Produção hidráulica  (Pn) MWh| 11367 12503 Basta LUSO ss as suga 52,3 84,9 
Produção térmica (Po) MWh 0 0 Deanaiã 19.8 99.9 
Produção total (PT) MWh/| 11467 | 12508 Pé bos is ao rm : 
Trocas com Export. (Ex) MWh Ú 0 OVOS e vv o» co qc. DX. 12,1 (3) 96,2 
Espanha | Import. (1)  MWh O O Total. . .| 818,5 18,5 
Total Pr -+(|-Ex) MWh| 11867 | 12503 , | 
Prod, para cons. perm. (Pep) MWh| 9458 10610 
Prod, para cons, não perm. (Pesp) MWh! 1909 1893 
ú — | Potênciamás, MW 6585 62 Notas : 
a Potência mir, MW 268,0 292 , ; p 
EE: Pr+ (FE | Teiliz. da ponta horas 17,8 18,1 (1) Coeficiente de enchimento em energia das albufeiras, 
À E a Sn | Of — Sib (2) Inclui 2,6 GWh armazenados em Vale do Rossim no 
ia Potência máx, MW 556 617 fim do mês 
É e p Potência min. MW 154 208 , 
So “Pp Utiliz. da ponta horas 17,0 17,2 (º) Inclui 2,7 GWh armazenados no açude do Poio no 
Factor de carga 0,71 0,12 fim do mês, 
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MAGUE 


CONSTRUÇÕES METALO-MECÂNICAS MAGUE, L.”* 
(S. E. MONIZ DA MAIA & VAZ GUEDES, L.”*) 
Alverca — Portugal 


Aspecto da maquinagem do ante-distribuidor de uma das 3 turbinas «Francis» de 120,000 CV destinadas 
à Central de Bemposta, construídas por Construções Metalomecânicas «MAGUE», 5. A. R, L., segundo 
projecto e licença dos Ateliers de Constructions Mécaniques de Vevey, S. A. 


“ 

PROJECTO E FABRICO DE 
PONTES ROLANTES 
GUINDASTES 
GRUAS. DERRICKS E GUINCHOS 
APARELHOS DE ELEVAÇÃO ESPECIAIS 
CONSTRUÇÕES METÁLICAS PESADAS 


TÉCNICA — XIII 


PRÉ-MONTAGEM 
DE 

UMA 
COMPORTA 

E DE 

UM 
RESERVATÓRIO 
NA 

OFICINA 

DE 
CONSTRUÇÕES 
METÁLICAS 


AÇOS VAZADOS AO CARBONO, MANGANÊS, CRÓMIO E INOXIDÁVEIS 
FERROS FUNDIDOS 

MECÂNICA MÉDIA E PESADA 

CONSTRUÇÕES METÁLICAS 


comportas 

turbinas 

condutas forçadas 

reservatórios sob pressão 

permutadores de calor — (Licença KELLOGG) 


REDUTORES DE VELOCIDADE 


DIVISÃO METALO-MECÂNICA 
COMPANHIA UNIAO FABRIL 


AV. 24 DE JULHO, 170 — LISBOA 2 — TELEFONES 670421 E 672061 — TELEGRAMAS: METALOFABRIL 


TÉCNICA — XIW 


C. D. U. 553.76 (469) 


Água fortemente radioactiva do Sul 


de Portugal Continental” 


A. HERCULANO DE CARVALHO 
C. MOURA PULIDO 


(Centro de Química — 1. S. T.—C. E. E. N. do Instituto de Alta Cultura) 


Numa nota publicada há cerca de 6 anos (1), 
um de nós tentou dar um panorama geral da 
riqueza radioactiva das nossas águas minero- 
-medicinais, tal como é conhecida através de 
trabalhos, próprios e alheios, trabalhos que, no 
departamento hidrológico do património nacio- 
nal são, felizmente, bastante numerosos, abran- 
gendo quase todas as nascentes de água mineral 
conhecidas no nosso território metropolitano. 

Uma das afirmações aí feitas era a de que 
todas as águas fortemente radioactivas (tomando 
arbitráriamente como limite inferior o teor de 
radão de 20 muC/1) se situavam a norte do Mon- 
dego. Verificava-se também que uma radioacti- 
vidade elevada das águas nem sempre estava 
ligada à existência, na região da nascente, de 
jazigos de minério radioactivo, isto é, que havia 
águas brotando fora das zonas uraníferas de 
Portugal Continental que apresentavam teor apre- 
ciável de radão (por exemplo, Estoril e Cucos). 
Pelo contrário, todas as águas minerais surdindo 
nessas zonas se mostram radioactivas. 

Os valores mais altos da concentração de ra- 
dão, conhecidos até então, situavam-se precisa- 
mente na região uranífera da Beira Alta, atin- 
gindo o nível de 50-60 mpC/l nas águas de 
Abrunhosa e S. Gemil, segundo as determina- 
ções de Egas Pinto Basto e A. Viana de Lemos (2). 


(1) A. Herculano de Carvalho. «Informação sobre ra- 
dioactividade das águas minerais portuguesas», Rev. 
Quím, Pura Apl., 36, 206, (1956). 

(2) «Determinações de Radioactividade em Águas Mi- 
nerais» — Rev. da Faculdade de Ciências da Universidade 
de Coimbra — Vol. VI, 1 e 3(1936-37). 


Apenas uma única nascente mostrou, aos mes- 
mos Autores, uma cota mais elevada e verda- 
deiramente excepcional: 179 mpC/1. 

Trata-se duma nascente, já assinalada e anali- 
sada por Lepierre, denominada «Urgeiriça» (1), 
situada a cerca de 300 m a sul do jazigo de urá- 
nio da mina do mesmo nome mas que, apesar 
de ter sido concessionada (1922) e duas vezes 
transmitida por alvará (1933 e 1948), nunca es- 
teve em exploração. 


Não eram pois conhecidas, ao sul do Mon- 
dego, nascentes de água mineral com tão eleva- 
das cotas de teor de radão. 

Recentemente, a convite da Câmara Municipal 
de Nisa, houve que proceder a análise duma 
água que nasce a cerca de 9 km daquela vila, e 
que, apesar de ser conhecida e utilizada, pelo 
menos desde os princípios do século passado, 
nunca fora examinada, sabendo-se apenas que 
era sulfúrea. 

As primeiras amostras para a determinação do 
radão indicaram uma ordem de grandeza muito 
alta do seu teor, o que aliás já se previa, dada a 
circunstância da nascente estar situada dentro 
duma zona em que os Serviços de Prospecção da 
Junta de Energia Nuclear tinham descoberto — 
facto inteiramente novo até há póuco, importan- 
tes jazidas de minérios de urânio. Supúnhamos, 
no entanto, encontrar um terço ou, quando 


(1) O valor indicado para a radioactividade refere-se 
à nascente «Urgeiriça II», mencionando os Autores cita- 
dos uma outra : «Urgeiriça I». 


* Comunicação apresentada à Classe de Ciências da Academia das Ciências de Lisboa em 12 de Julho de 1962, 
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muito, metade do valor que depois foi determi- 
nado e, por isso, as medidas não se fizeram em 
boas condições de precisão, dada a sensibilidade 
do dispositivo de medida. 

Assim, foi necessário voltar ao lugar da nas- 
cente, repetir as colheitas e, já em boas condi- 
ções, as medidas. 


A nascente, situada entre a sede do concelho 
e as povoações de Arez, Tolosa e Gáfete, é 
conhecida pelo nome de «Fedegosa» ou «Fada- 
gosa». Esta designação, justificada pelo cheiro 
das águas, é comum a mais três mananciais de 
água sulfúrea na mesma zona: a de Mação, no 
distrito de Santarém, e as de Gavião (l) e de 
Marvão (também chamada Maria Viegas ou Pe- 
reiro) ambas no distrito de Portalegre. Os locais 
das quatro Fedegosas podem situar-se dentro 
dum rectângulo, de lados paralelos às direcções 
oeste-leste e norte-sul com as dimensões de 55 
e 15 km, respectivamente. 

Em «Le Portugal Hydrologique et Climatique» 
encontra-se uma análise (Lepierre-1932) da Fada- 
gosa de Gavião. A água da Fadagosa de Mação foi 
resumidamente analisada pelo Dr. Mastbaum e, 
posteriormente, por um de nós; quanto à água 
de Maria Viegas, apenas se conhece uma análise 
muito sumária, de 1891, devida a Santos e Silva. 

Como notámos, nada se sabia até agora sobre 
a composição da água da Fadagosa de Nisa ou 
de Arez. 

A análise mostra tratar-se duma água fria, 
do tipo bem caracterizado das sulfúreas-alcalinas 
e portanto hipo-salina, muito fluoretada (10,8 mg/1 
de F ), com predominância muito marcada do 
catião sódio e do anião bicarbonato. 

Tem reacção alcalina (pH = 7,92) e é redutora 
(e red =— 0,12 V; rH: = 11,8) (2). A não ser 
por um teor um pouco mais baixo de sílica, não 
se distingue das suas congéneres do Centro e 
Norte do país. 

Mas o que verdadeiramente a individualiza é 
a sua radioactividade. Esta característica, em 
absoluto, e a circunstância dela se apresentar, 
pela primeira vez até hoje, numa água do sul 


(1) Hoje inutilizada por inundação após a construção 
da barragem de Belver. 

(2) Determinações de J. Duarte de Almeida, do Insti- 
tuto de Hidrologia de Lisboa. 
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de Portugal Continental, justifica a apresentação 
da presente comunicação. 

Há pouco que dizer quanto aos métodos de 
análise que, quanto às determinações do radão 
e do rádio, são os habituais no Laboratório do 
I. S. T., onde trabalham alguns bolseiros do 
Centro de Química Nuclear do 1. A. €. 

Na determinação do urânio, usou-se o método 
muito sensível de separação, por extracção com 
benzeno e dosagem espectrofotométrica do com- 
plexo com «rodamina B», método estudado e 
desenvolvido neste Centro de Estudos (1). 

Os resultados condensam-se na tabela a seguir: 


Por Por f 
| 1 litro |1 gderesiduo Unidades 
DS in 
Radão (Rn) | 84,2 | — | mel 
Rádio (Ra) 1 3 x 10-1º g 
Urânio (UsOs) | 135,6 410 Kg | 
| 


Vê-se assim que não só o teor de radão é ex- 
cepcional (o mais alto de todas as nossas águas 
em exploração) mas ainda que: 


a) a concentração de rádio se situa para lá 
do limite considerado normal (107 4!); 
b) o teor de urânio é o mais alto de que 
tenho conhecimento em águas naturais 
portuguesas, cerca de 100 vezes o valor 
que usualmente se encontra. Sob este 
aspecto, observa-se que o resíduo da 
água de Arez poderia considerar-se já 
como um minério de urânio pobre (cerca 


de 0,05 “/, de UsO)). 


Outras considerações merecem ser feitas, opor- 
tunamente, sobre certos aspectos químicos da 
coexistência destas características numa água 
natural. 

Por agora, queríamos apenas pôr em evidência 
a raridade do facto, propondo aos médicos e aos 
geólogos um tema de meditação sobre as suas 
implicações nos planos que são do âmbito das 
suas especialidades. 


(1) Fraústo da Silva e M. Legrand de Moura — Novo 
método espectrofotométrico para dosagem do Urânio — 
Técnica, 26, 258, (1958). 


C. D. U. 536.24.004 


Condução do calor em sólidos em regime variável 


(Soluções pelo método de separação de variáveis) 


PELO ENG.º MEC. JOSÉ JOAQUIM DELGADO DOMINGOS 


Assistente do 1.5. T. 
Bolseiro da OTAN 


À importância da condução do calor em sólidos como um dos modos fundamentais 
da transmissão do calor, junta-se a importância crescente do estudo das tensões em sólidos 
provenientes de gradientes térmicos. Se sob aquele aspecto as soluções analíticas têm interesse 
na engenharia, sob este essa importância acentua-se por ser um dado primário do problema. 
Julga-se por isso de interesse uma sucinta tentativa de correcta definição do problema da 
condução do calor em sólidos do ponto de vista macroscópico por forma a situar os modelos 
usualmente utilizados na sclução analítica. Como exemplificativo apresenta-se um dos 
métodos da análise mais conhecidos e utilizados na obtenção de soluções — o método de sepa- 
ração de variáveis. Pensa-se de interesse a apresentação da bibliografia mais importante, bem 
como a resumida indicação das vantagens e limitações da utilização do método nas aplicações 


O — Introdução 


0.1 — O estudo da condução do calor faz-se 
com toda a generalidade em Termodinâmica Irre- 
versível na qual, sob o ponto de vista macroscó- 
pico, se postula a validade das funções termoestá- 
ticas relativamente aos valores instantâneos das 
variáveis e se admite acoplamento entre fluxos e 
afinidades diferentes. Se o sistema se não encon- 
tra muito afastado do equilíbrio termoestático 
aquele acoplamento pode exprimir-se por rela- 
ções lineares sob a forma duma matriz de inter- 
ligação, que é simétrica no caso da inexistência 
dum campo magnético exterior, e os seus elemen- 
tos são no caso geral função de todos os poten- 
ciais. Além disso, como a soma em cada ponto (1) 
dos produtos de cada fluxo pela respectiva afini- 
dade é igual a produção da entropia nesse ponto 
por unidade de tempo devido à irreversibilidade 
aquela soma é positiva e a matriz é defenida 
positiva. 

Apesar do seu interesse, principalmente quando 
se trata o estudo de sistemas com acoplamento 
de fluxos, não é corrente a utilização da Termo- 


(1) Entendido no sentido físico, i. e. num elemento 
de volume macroscôpicamente muito pequeno. 


dinâmica Irreversível no domínio da engenharia 
e não se encontra nos tratados recentes da lite- 
ratura ocidental no domínio da transmissão do 
calor e de massa (2), 


0.2 — Não se considerando a existência de 
acoplamentos de fluxos, a forma clássica de abor- 
dar a condução do calor é a partir da equação 
de Fourier, a qual é, relativamente à condução 
do calor, o equivalente da lei de Ohm em elec- 
tricidade. 

No contesto da Termodinâmica Irreversível 
significa que apenas são diferentes de zero os 
elementos da diagonal principal da matriz, ou 
seja, que os vários fluxos só dependem das res- 
pectivas afinidades e além disso que cada ele- 
mento da diagonal é apenas função do respec- 
tivo potencial. Tal permite que os vários fluxos 
e respectivas afinidades sejam estudadas inde- 


(2) A Termodinâmica Irreversível que é de cria- 
ção essencialmente europeia, com os trabalhos de 
Onsager, Prigogine e De Groot, principalmente, é utili- 
zada neste domínio pelos autores orientais, em parti- 
cular Lykov e Mikhaylov na sua Teoria da Transmis- 
são da Energia e Massa em que a aplicam extensa- 
mente no estudo da secagem. 
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pendentemente dos restantes, e estudar portanto 
isoladamente a condução do calor (8). 

O ponto de vista clássico corresponde a esta 
simplificação. No entanto, embora isolando a 
condução do calor e mantendo sempre a ausên- 
cia de acoplamento de fluxos, a teoria clássica 
têm-os em conta em certos casos introduzindo 
no sólido fontes ou poços de calor, ou artifícios 
equivalentes, pelas quais considera a produção 
da entropia pela irreversibilidade a elas associada 
(v. nota n.º 7) (3), 

Esta simplificação justifica-se quase sempre 
nas aplicações porque reduz-se a complexidade 
do problema o qual do ponto de vista matemá- 
tico passa a exprimir-se, no caso geral, por um 
sistema de equações às derivadas parciais inde- 
pendentes em vez de o ser por um sistema de 
equações simultâneas mais complexas. Obvia- 
mente, considerando apenas a condução, o pro- 
blema reduz-se a uma só equação às derivadas 
parciais. 


0.3 — Além das simplificações anteriores é 
usual considerar um sólido homogéneo e isó- 
tropo e admitir a constância relativamente à 
temperatura da constante de proporcionalidade 
da equação de Fourier. Nessas condições se obtém 
a conhecida equação às derivadas parciais do tipo 
parabólico conhecida por equação do calor: 


1 OT 
x dt 


a qual, como se verifica, corresponde a uma sim- 
plificação do modelo físico mais geral que se 
esboçou. 


0.4 — A simplificação que conduz à equação 
do calor é inevitável em face da complexidade 
analítica que resultaria da consideração do mo- 
delo mais geral que a Termodinâmica Irreversi- 
vel oferece, 

Além daquela complexidade, outra objecção se 
levantaria nas aplicações e que seria resultante 
da imprecisão com que se conhecem as caracte- 
rísticas térmicas dos materiais, ou os constrangi- 


(3) De igual modo terá em conta a existência dos 
potenciais nos elementos da matriz defenindo estas o 
valor constante do potencial que seja necessário con- 
siderar. 
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mentos reais que se exercem no sólido, os quais 
são característicos do problema a resolver. 

No domínio corrente das aplicações é portanto 
geralmente admissível a consideração de modelos 
simplificados e é de alguns deles que nos iremos 
ocupar tendo esta curta introdução apenas a fi- 
nalidade de colocar o problema dentro dum con- 
texto mais geral que permita melhor situar os 
casos que usualmente se tratam. 


1 — Formulação do modelo 


1.0 — Simbologia 


T, T* — Temperatura — (a origem da escala é 
arbitrária) 

Ti e To — Temperatura inicial do sólido e do 
meio exterior, respectivamente — constantes 

t — tempo 

U — energia interna 

Q — quantidade de calor (por unidade de vo- 
lume e tempo — fonte — nalguns casos) 

c — calor especifico 


— . 


q — fluxo térmico 
o — massa específica 


K — condutividade térmica 
K 


« — difusividade térmica == Fi 
k 


h — coeficiente de transmissão de superfície 


E A at 
Fi — Número de Fourier - cu 


B; — Número de Biot — h l 
H = h 
K 


P — ponto genérico no sólido 

x, Y, Z — coordenadas cartezianas 

r — raio (coordenadas esféricas ou cilíndricas) 
R — raio da superfície exterior 

D — domínio em que existe o sólido 

B — contorno de D 

Bc— parte de B comum a dois domínios 


n — normal exterior a B 

V — volume 

S — superfície 

| — metade da espessura da placa 

q”, É, y — valores próprios 

v — função própria 

WY — função do tempo 

F (P) — temperatura no sólido para t == 


[] — bibliografia 
() — notas à margem 


— Formulação 


1.1 — No estudo da condução do ealor em só- 
lidos admitiremos que: (4) 


a) O sólido em estudo está contido numa 
região D limitada por um contorno B; 

b) O calor é a única forma de energia per- 
mutada entre as várias regiões do sólido 
ou deste com o exterior através do con- 
torno (5) ; 

c) Em cada elemento de volume do sólido 
a variação de energia interna se exprime 
por 


dU=c.p.dT (6) 


d) Existirão geralmente no sólido fontes 
(ou poços) de calor de intensidade Q (70) 
por unidade de volume (7); 


(4) Supõe-se implícito o tratamento macroscópico 
e o sólido contínuo do ponto de vista matemático. 
Tal implica a consideração de elementos de volume 
macroscôpicamente muito pequenos mas microscôpi- 
camente muito grandes. Prâáticamente a distância mí- 
nima que se considera é superior às distâncias interató- 
micas. 

(5) O que implica que o sólido se não deforma ou 
que o trabalho a este associado é nulo — em rigor; 
em geral considera-se que é desprezável, 

(6) Não se considerando outras formas de energia 
que não calor, e tendo em conta b) c identifica-se com 
cy — calor específico a deformação constante. Como 
porém em geral embora desprezando a deformação 
esta não é nula, c não fica bem defenido. Este ponto 
encontra-se omisso em quase todas a tratados sobre 
a condução do calor, nos quais quando se identifica c 
é com cp. Tal podé fazer-se supondo nulo o trabalho 
interno de deformação, e que o sólido troca trabalho 
com um meio exterior (através dé B) a pressão cons- 
tante embora sem alterar prâticamente a sua forma 
geométrica. 

Prefere-se a formulação adoptada porque com- 
preendendo as anteriores permite a melhor identi- 
ficação consoante os constrangimentos reais — conti- 
nuando ôbviamente a considerar sempre apenas o caso 
em que não há alteração apreciável das formas geo- 
métricas. (Aliás num sólido, e dentro destas condições 
cp e cy são práticamente iguais numéricamente). 

(7) A consideração de fontes de calor — fontes de 
entropia - é geralmente como já se indicou a forma 


e) Não se verificam mudanças de fase. 

f-) Entre a temperatura e o fluxo de calor 
num ponto verifica-se a relação de Fou- 
rier (8), 


q K grad T 
sendo K > 0 


1.2 — Dentro das hipóteses anteriores se con- 
siderarmos um sólido de volume V, limitado 
pela superfície S, a conservação da energia 
implica que: 


O dV= [= dV 
e dt 


— [a.nds+ 
s V v 


ou seja 


- fava av + [odv= [Say 
V “V v 


———————— e ————— 


simples de desprezando a existência de outras trocas 
de energia (não calor) se ter em conta a irreversibi- 
lidade dessa permuta. Por exemplo, não se terá em 
conta num condutor percorrido por uma corrente eléc- 
trica o efeito do acoplamento deste fluxo, ou de exis- 
tência do respectivo potencial mas dar-se-á conta da 
irreversibilidade que se traduz pelo efeito de Joule 
introduzindo este sob a forma duma fonte de calor. 

Também se consideram por vezes fontes referidas 
a unidade de superfície ou de comprimento. Embora 
sem realidade física, tal como oanterior, justifica-se por 
simplificação do modelo físico ou por artifício analítico, 

(8) Implicitamente admite-se a isotropia do sólido. 
Se o sólido não é isótropo K é um tensor de segunda 
ordem simétrico. Neste caso cada componente do fluxo 
térmico segundo uma direcção é função linear das 
componentes do gradiente, os coeficientes sãotermos 
da matriz de ligação entre fluxos e afinidades, (neste 
caso referentes a fluxos que embora diferentes são 
da mesma natureza) — V. [5], [6]. 

A anisotropia raramente se considera nas aplicações. 
Ve. [x], [3], (4). 

— A relacão de Fourier, implica uma transmissão 
instantânea da perturbação. Tal não se verifica, mas 
admite se porque a velocidade de propagação é 
muito grande. Tendo-a em conta, introduz-se em 


= 
dq 
ot 
sendo 3 uma constante de tempo. Com esta relação 
obter-se-ia não a equação do calor mas a equação dos 
telegrafistas, a qual só pode vir a ter interesse quando 
se considera a condução do calor em gases rare- 
feitos p. ex. 


primeira apoximação no primeiro membro à a 
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donde passando ao limite, será 


E diva + Q — 9U 
dt 
e por c) e f) 
9T 
div (— K grad T) + do E 
t 
ou, usando o operador yv 
11 v(KyT)+Q=ee 
t 


se o sólido é homogéneo e K independente de T, 


1.1 pode escrever-se fazendo E = « — difusi- 


vidade térmica 


12 a A Q 
a dt ão 
ese Q=0 
1.3 vºT E ET 
qa dt 


e se o regime é estacionário 
1.4 VET = 


A equação 1.3 é a conhecida equação do calor. 
1.4 é a equação de Laplace. Se o regime é esta- 
cionário mas Q + 0 


1,5 4T-—= — Q 
a RE 


que é a equação de Poisson. 


1.3 — Na sua forma geral — 1.1, a equação não 
é linear nem homogénea e a solução analítica é 
difícil de obter (9), 

de K é apenas função da temperatura pode 
obter-se uma equação formalmente análoga a 
1.2 introduzindo uma nova escala de tempe- 
raturas por 


à PE 
TT) = - 
(T) e). KM AT 


com o que 1.2 se pode escrever 


1.2 Tres E ZE E: 


(9) Para o caso geral conhecem-se muito poucas 
soluções analíticas mesmo para formas simples de 
K(P, T) e Q(P). 
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Embora a consideração de K = K(T) tenha 
importância sempre que se considera um inter- 
valo apreciável de temperaturas, K admite-se 
geralmente constante. A variação de K com T 
admite-se, quando se tem em conta, quase sem- 
pre linear, mas as aplicações são em geral difi- 
cultadas porque mesmo nesta forma simples os 
coeficientes de variação apenas se conhecem para 
pouco materiais (10), 


1.4 — As soluções analíticas que se conhecem 
em maior número referem-se à equação 1.3. 
Dado o seu interesse nas aplicações é dela que 
em particular nos ocuparemos enquanto conside- 
rados os métodos analíticos. 


2 — Condições iniciais e no contorno 


Os problemas em condução do calor sinteti- 
zam-se em determinar : 


a) a função que satisfaz a equação às deriva- 
das parciais dada, que 

b) Verifique certas condições no contorno 
(superfície no sólido) e que 

c) Assuma num dado constante (tomado para 
inicial geralmente) certo valor no domínio. 


As condições a), b), c) garantem a unicidade 
da solução (11) como o exige o sentido físico do 
problema [8]. 

Dentro desta formulação geral cada par de 
condições (iniciais e no contorno) representa um 
problema particular, o que é uma caracteristica 
dos problemas físicos que se formulam matema- 
ticamente deste modo. Existe assim um número 
infinito de problemas de que se não conhece 
método analítico de obter a solução. Por este 
motivo e em ordem a simplificar o tratamento 
analítico apenas se considera usualmente um 
número restrito de condições no contorno, e se 
admite que T em D satisfaz as condições de 
Dirichelet. 


(10) Será objecto de futura publicação o método 
que seguimos no N. E. E. M. para a sua determinação. 

(11) À unicidade de solução demonstra-se matemã- 
E , a, EL e a. limitados 

dm dx dy dz 

em Be D respectivamente e T contínuo em t(— t==0) 
[1]. No caso geral o problema de unicidade não tem 
interesse nas aplicações, em que, quando necessário é 
testada em cada caso particular tratado. 


ticamente para 


2.1 — Condições no contorno 


As condições que é normal considerar no con- 
torno (B) são: 


2.1.1 — Temperatura em B função imposta de 
ter 


2.1 Tp = T (P, t) B 


2.1.2 — Fluxo de calor imposto através de B 


Kk—— | = fá 
( dn ). a ), 


ou, como caso particular, e correspondendo ao 
caso de isolamento térmico perfeito: 


2.3 (x 


2.1.3. — Fluxo de calor através de B propor- 
cional à diferença de temperatura entre Beo 
meio exterior suposto a To uniforme. 


2.4 (x = lh LTo—T (2,0 ) 


n 


2.2 


B 
eh > 0 (12) e geralmente h = h (P, t), 
h será designado por coeficiente de transmis- 
são de superfície (se h = 0 2.4 transforma-se 
o T 
em 2.3; se h=-0202, como K a é finito por 


n 
hipótese, deverá ser To = T, e 2.4 transforma-se 


em 2.1). 


2.1.4 — Dois sólidos 1 e 2 em contacto ao longo 
de Bo sem resistência térmica de contacto: 


Ta (P, De = Ta (P,t) 
d 13) 
2.5 (xa ) E (x: eg ( 
on Be dn Bc 
Se existe uma resistência de contacto R a 
1.2 equação 2.5 é substituída por uma condição 


do tipo 2.4 com h = a e To = Ts (Pe,t)p. 


2.1.5 — Apesar de formalmente simplificado o 
tratamento analítico que se conhece impõe ge- 


(12) Esta condição resulta simplesmente do 2.º prin- 
cípio da termodinâmica pois K >o. 


e 


— —» 
(13) Será n=n, ou n=n; a normal comum aos 
dois sólidos sobre Bc — sendo óúbviamente exterior 
em relação a um e interior relativamente ao outro. 


ralmente restrições importantes nas condições 
anteriores. Assim, h e as restantes funções são 
consideradas constantes na superfície (B) ou 
pelo menos em domínios finitos em que aquela 
se subdivide — quanto à dependência do ponto. 
Quanto à dependência no tempo o problema 
trata-se recorrendo ao teorema de Duhamel que 
reduz o problema a um de condições no con- 
torno independentes do tempo (14), 


3 — Métodos analíticos em regime variável 


3.0 — Como se indicou a forma mais geral que 
usualmente se considera da equação do calor 
é 1.2. 

Em face do anterior o problema formula-se 
pois, matematicamente, do seguinte modo: 


Determinar T=T(P,t) para P em B+D 
e t20. 


Tal que 
ap 1ºTl O > 
ve T E = para P em Det >0. 


Sujeita a uma das condições 2.1, 2.2, 2.3, 2.4 
ou 2.5, ou sua combinação (15) para P em B 
e tal que 


T(P,0)=-F(P) para Pem De t=0 


Não se tratando da consideração exaustiva 
dos métodos matemáticos de tratamento do pro- 
blema — que é um dos capítulos mais importan- 
tes da análise — abordaremos apenas, alguns dos 
mais úteis do ponto de vista do engenheiro, que 
no presente trabalho se reduzirão ao método de 
separação de variáveis e proximamente os das 
transformadas integrais, como exemplificativas 
de métodos analíticos. 


3.1 — Método de separação de variáveis 


3.1.0 — Consideremos as seguintes simplifica- 
ções relativamente ao modelo mais geral pro- 
posto : 


a) k,e, c são constantes (e .'. 2) 


(14) V. [17] para a demonstração clássica. Será por 
nós apresentado quando tratarmos da aplicação do cál- 
culo operacional, onde este teorema se identifica com 
o integral de Faltung, ou da convulsão [9]. 

(15) Obviamente B poderá subdividir-se em sub- 
contornos Ba, Bb, Bc...Bj, e B=IBij. 

Sendo em cada B; aplicável uma condição. 
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bh) o domínio D+B é finito e regular 
o) F(P) satisfaz em D às condições de Diri- 
chelet 
d) T(P,t) é de classe C(2) em (D+ B) 
para t > 0) 
e) sobre o contorno (B) verificam-se condições 
do tipo 2.4 independentes de t (16) 
f) Q(P) = 0 (16) 
e admitamos, o que é 
que 
g) Ti; (P, t) satisfazendo a equação do calor e 
as condições no contorno se pode exprimir 
em D pelo produto duma função só de P 
e outra só de t não idênticamente nulas (17) 
para t 20 
Ti (Pt) = qa (P) 4 (t) 
h) o sistema de funções próprias 7; se orto- 
gonal, ou ortogonalizável constitui um 
sistema completo [8, 9]. 


é característico do método, 


3.1.0.1 — Tomemos como origem da escala de 
E já 
temperaturas To e façamos H = E -(>0). 


Pretendemos determinar T (P,t) tal que: 


1.3 vT— E a = 0 PemDet>D>o 
a dt 
7 
2.4 E ce PemBet>oO 
n 


31 T(P,0) = F(P) PemBet=-0 


Como por hipótese: 


3.2 T;i= o(P) 'Y(t), satisfaz 1.3 e 2.4 
obtém-se, substituindo em 1.3, que deve ser 


vel) 1 Tt) 


3.3 ee 
2 (P) a Y(t 


(2 > 0) 


o 1º membro é só função de Pe o segundo só de t 
— são pois iguais a uma constante, que designa- 
remos por 2. 

Será portanto 


34 va(P)-MY(P)=0 PemD 
W ; 

3.5 É — haVf=0 t>Do 
dt 


(16) V. capítulo 3.2 condições mais gerais. 

(17) Considerar alguma ou ambas idênticamente 
nulas implicará admitir-se um equilíbrio termoestático 
indefenidamente. 
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Substituindo 3.2 em 2.4” e atendendo a que 
esta derivação só envolve coordenadas espaciais, 
vem 


) 
vo [Ens 1|=o paraPem Bet >o0 
e como W(t) não é idênticamente nulo, por 
hipótese, deverá ser 


GR 
3.6 


= Elgg mes 


para P em B 
dn 


3.1.0.1.1 — Considere-se agora o teorema de Ostro- 
gradsky : 


f v.adv=[ (aa) ds 
v “ g 


sejam fi e fz duas funções escalares do ponto e 
faça-se a == fi y f: obtém-se, substituindo aci- 
ma, a identidade 


37 [vitor th ii die =[ hitds 


Façamos agora fi=f)=-%, obtém-se da 
anterior (18) fazendo V=D a S=B 


—— — jp d 
favicav+ INCZORITO dv=[ a ds 
D JD B n 


os 
substituindo v2v de 34 e E de 3.6 resulta, 


atendendo a que 2; é constante e 7; não idênti- 
camente nulo: 


8 fis: (v 4). (v 3) dv + Í, H (Ga) 
fo 


como H>>0 por defenição, a fracção é não 
negativa e é portanto 


3.9 M4A<O 
E o sinal de igualdade só poderá verificar-se 


do 
se H=0 ou —!=0 e V9WÔ=o. 


dn 
Esta última condição significa que 7; == const 
e portanto que a temperatura se mantém indefe- 
nidamente constante, pois Y é não idênticamente 


(18) Por d) e 2) e pela conclusão a que se chega 
quanto a 4, « satisfaz as condições para ser válida 
a identidade. 
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nulo e é limitado como se verá. A condição 
0 o : 

H =0 ou ra = O significaria um sistema iso- 
n 

lado donde a conclusão física trivial de que o 

sistema em equilíbrio Termoestático permanece- 

ria indetenidamente nesse estado porque isolado. 

Para acentuarmos o carácter não positivo 
de 4; faremos 


3.10 


o 


h=— q 


considere-se agora 3.5; integrando e atendendo 

a 3.9 obtém-se 
3:11 

sendo c uma constante arbitrária (HF 0) 

de 3.11 conclui-se que '» é limitada (pois q) > 0) 

como se afirmou (19), 


= 2 at 
Wc q; 


3.1.0.1.2 — os q; são como se sabe os va- 
lores próprios da equação de Helmoltz no pro- 
blema defenido pelas condições no contorno 3,6. 
A cada um corresponde uma tunção própria 
cy (20). Os q%, para que %; seja função própria 
não são quaisquer, podendo assumir apenas um 
número, discreto de valores (21), 

Verifiquemos agora que o sistema de funções 
próprias 2; . .. 7ycorrespondentesa q;.... q; são 
ortogonais. 

De 3.7 permutando fi; e f: e subtraindo as 
duas expressões obtém-se a identidade de Green: 


d 
fiva—tvtt) dv = f(a E — 
D B 


n 
mods 


dn 


(19) Tal corresponde à conclusão física, evidente, 
de que no sistema isolado constituído pelo sólido 
-+- meio exterior a entropia não pode diminuir. A forma 
clássica de justificar a afirmação à priori de quez;<,o 
deriva directamente desta conclusão, que se exprime 
dizendo que por razões de ordem física T não pode 
crescer indefenidamente. 

Chama-se contudo a atenção para o facto desta con- 
clusão ser obtida sem necessidade de introduzir outras 
hipóteses além das referidas na formulação do modelo. 
Na forma usual de proceder, existe redundância nos 
postulados, como se vê. 

O ser »j <, o corresponde a afirmar que o sistema 
é puramente dissipativo — perturbado tende para o 
equilíbrio por relaxação pura. 

(20) De notar que a um q; pode corresponder mais 
do que uma função própria — o caso aflora-se adiante. 

(21) Não se faz a demonstração [8 - 10] — supõe-se 
na formulação contida em afirmar que «j é função 
própria. 


—— Pa 


fazendo 7; = fi e 7; = f> resulta atendendo nova- 
mente a 3.4 e 3.6 


[in adn—s. ; n) dyv= 
o D 


“e 


fes ] (—Hgg+hH ga) ds 
B 


e por serem q; e q; constantes 


312 (q— 95) ) 4 dv=o0 
D 


donde, se q; * q; 


fa 2 dv=0 
D 


o que exprime a ortogonalidade das funções 
próprias no domínio D (22), 


3.13 (15=1) 


3.1.0.1.3 — Demonstradas as características 
particulares de v e 'Y para que T; = u(P) W(t) 
seja solução de 1.3 e 2.4” resta encontrar a fun- 
ção que satisfazendo também 3.1 é a solução do 
problema. 

Ora, dada a linearidade das equações 3.1 e 2.4 


to Por 2a 
se T; for solução Ts T = C; x (P) e q 2 


também é solução de 1.3 e 2.4' se a série for 
convergente. 
Para que T(P,t) seja solução do problema é 
necessário que T(P,0)=2C;z (P) =F(P). 
| 


Mas o sistema das funções próprias é ortogo- 
nal (ou ortogonalizável) e por hipótese é com- 
pleto. Segue-se que F(P) pode ser representado 
em D pelo desenvolvimento em série dos %;, 
para o que basta determinar os coeficientes C,. 
Estes determinam-se atendendo às relações de 
ortogonalidade : 


[EM auav=i | Genmudy=6 [st dy 
D k o ” D 


(22) Como se indicou, a um valor próprio podem 
corresponder funções próprias diferentes — por ex. 
qi e di — e neste caso 3.12 não demonstra a ortogona- 
lidade pois o parênteses é nulo. Mas os « e &j são li- 
nearmente independentes, como é óbvio, o que per- 
mite construir a partir delas por combinação linear um 
sistema ortogonal. Dada a linearidade das equações o 
sistema é equivalente. Por se afastar demasiado do 
âmbito do trabalho não se demonstrará aqui este as- 
pecto particular — V. p. ex. [8). 
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| F (P) q; dv 
3.14 C; = e Lo 


] Dj dv 
D 


e a solução do problema será 


os —g2a 
315 T(IL)= 3 Gau(Pe 
|=1 


3.1.0.2 — Exemplo analítico. 


Considere-se o problema clássico duma placa 
inflnita segundo y e z, de espessura 21 (—l<x<]1), 
inicialmente à temperatura T (x, 0) == EF (x), ins- 
tante em que se mergulha num meio a tempe- 
ratura To, constante e uniforme. O coeficiente de 
superfície h é igual para ambas as faces e cons- 
tante, e a placa supõe-se homogênea e isótropa, 
e de propriedades térmicas (Kec) independentes 
da temperatura. Pretende saber-se a evolução 
das temperaturas no interior da placa ao longo 
do tempo. 

Como a temperatura inicial é só função de x, 
e como a sua extensão infinita segundo y e z 
equivale a um isolamento segundo planos perpen- 
diculares a essas direcções, a temperatura deverá 
continuar a ser apenas (por razões de simetria) 
função só de x. O problema trata-se matemáti- 
camente a uma dimensão, sendo o domínio D 
—I<x<<leo contorno x =+1. 

A sua formulação é, tomando To; como origem 


h 


da escala da temperatura e pondo H = Te 


Determinar T ==T (x,t) tal que 


à 
3.16 E çã E) para —-|I<x<l 
dx? a dt 
e + >» 0 
e que 
3.17 E HT O ql (sô 
Ox /+1 
sas (05) = .HT(=1,t) x==| t>0 
OX /—1 
Sendo 
349 T(ix,0)=E(o) -I<qx<1l e t=0 


A solução obtém-se pelo método indicado que 
se aplica agora directamente. 
A equação de Helmoltz correspondente obtém- 
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-se escrevendo em coordenadas cartesianas e a 
uma dimensão 


3.20 aT e qu 1 ==) 
dx 

Trata-se neste caso duma equação deferencial 
ordinária, o que permite individualizar imedia- 
tamente as possíveis funções próprias que esta- 
rão óbviamente contidas no integral geral (23), 


3.21 FK(x)=Acosqx+B.senqx 


de que podem ser os seguintes casos particulares 
(exceptuando o caso de 3.21 já o ser) com 
| Ai + |Bil & 0). 


a) Ao BED 
b) A= 0 B50 


Verifiquemos sucessivamente a) e b). 


3.1.0.2.1L — Para que f=- A cos q x (Ai 0) 
seja função própria deve verificar 3.17 e 3.18 ou 
seja, deve ser simultaneamente : 


Aqsenql=AHCcosql 
Agsen(—-ql)=— A Hcos (— ql) 


Como se verifica as duas equações conduzem 
à mesma condição (Pi -& 0). 
qsenqil=-Hcos ql 
que se costuma escrever 
3.22 (gltag (q]) = HI] 
igualdade que só é satisfeita para os valores dis- 
cretos do parâmetro q; (i=-=1...n) que são 
raizes da equação transcendente 3.22. Estas raízes 
podem determinar-se graficamente verificando-se 
com facilidade que são em número infinito (24) 
e reais 
Considera-se agora o caso b): 
= B senqx iBj0) 
Analogamente se conclui que os q não são 
quaisquer mas apenas as raízes da equação trans- 
cendente (24) 


3.23 (qDcotg(gq)+ Hl=o 


(23) Se se tratasse do caso geral da equação às de- 
rivadas parciais a solução geral só raramente teria 
interesse prático dadas as suas características muito 
genéricas. Começa-se nesse caso, usualmente, por pro- 
curar funções que satisfaçam as condições no con- 
torno. V. [5] p. ex. 

(24) Em [1] encontram-se tabeladas as primeiras 
seis raízes. Em [2], [4], [13] também podem encon- 
trar-se tabelas, mas mais incompletas. 


que são reais e em número infinito qr (K =1....). 
Verifica-se também que sendo qi= q: 3.21 
não poderá ser função própria. 
Ponhamos 


para distinguirmos os valores próprios corres- 
pondentes a cada caso. O sistema de funções 
próprias será pois 


cos &, sen dy, cos És, sen ds... 


que já sabemos ser ortogonal. É além disso com- 
pleto para o que basta notar que se identifica 
com os termos duma série de Fourier. 

Se F(x) satisfaz as condições de Dirichelet, 
no intervalo —| <x<1 é nesse intervalo 


F (x) =s (A; cos Bix + B; sen 0; x) 
-=1 


Obtendo se por 3.14 os coeficientes: 


adel 
| - FG) cos (Ex) dx 
3.24 Aj= Ns 
] cos? (B;x) dx 
+! 
f F (x) sen (2x) 
a fe E ecenççeo 


.. i 
sen* (0;x) dx 


donde a solução do problema 
oo mi 2 at 
3.26 T(x, ps [4 cos (ix) e ! + 
==] 
" — Bjat 
+ B;sen (:;x) e 


Havendo apenas que juntar T, ao segundo 
membro para que T se exprima na escala de 
temperaturas onde o meio exterior está a To. 


3.1.0.2.2 — Consideremos o caso particular da 
temperatura inicial ser uniforme e igual a Ti 
relativamente à origem da escala em To. 


Será então F(x) = T;= const=-T(x,0) 


como x ==0 coincide com o centro da placa 


segue-se que F (x) é par donde B; = O (25), 
E obtém-se para A; 


p+1 
Ti] cos (f; x) dx 
3.27 A=D2—— — = 
f A, cos” (3; x) dx 
— 7 2sen (Bl) o 
sen (2; 1) cos (É; 1) + il 
donde 
3.28 T(x,t) = 
a o 
= Ti E doses! é II cos (14) 
i=1 sen cos 7 + 7i va 
tendo-se feito para simplicidade de escrita 
n = El 


se considerarmos a escala habitual de tempera- 
turas, teremos, de 3.28 e por defenição de T e Ti: 
3.29 ES Pi 

Ti—T, 


at 
— 4% —— 


|? 
x 
cos (1 8) 


A expressão 3.29 é adimensional o que se presta 
particularmente bem à representação grática das 
soluções o que se faz usualmente tomando os 
dois números adimensionais 


o . 
, sen 


i=l sen) cos 7; + 7 


B;=H1l = e — número de Biot 


at 


—— — número de Fourier 


Fu 
É 


um como parâmetro o outro como variável (26), 
e considerando x ==0 oux ==1,e por vezes gráficos 


à x 
para a conversão para qualquer valor de — 


(29) Deve notar-se relativamente ao problema de 
E (x) não constante que os valores próprios são função 
da «situação» da placa relativamente ao eixo das coor- 
denadas. Se a tivessemos defenido por exemplo no 
domínio o < x <,21 obtinha-se que B="o0 e T(x,t) 
exprimia-se apenas por uma série de cosenos. Prefe- 
riu-se o domínio — Il <,x <,l para que a solução 
exemplificasse o mais amplamente possível o método. 

(26) Tais representações, efectuadas também para o 
caso do cilindro, e esfera foram feitas entre outros por 
Gurney-Lurie [14], Gróber [4], Hottel [14], Pachkis [2], 
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[15, 14... | Como fácilmente se verifica a série 
converge rápidamente com o número de Fourier 
e como os 7); são crescentes podem geralmente 
usar-se apenas os primeiros termos para números 
de Fourier pequenos. 

Para valores muito pequenos do número de 
Fourier, o que interessa por exemplo no estudo 
do choque térmico, tratamentos térmicos, etc., a 
cérie converge lentamente e é pouco prática a sua 
utilização. Deve então usar-se outra forma de 
solução que é equivalente nas rapidamente con- 
vergente para F; pequeno e que se obterá quando 
abordarmos o cálculo operacional. 

Deve notar-se, como já se acentuou, que o caso 
de a condição no contorno ser uma temperatura 
dada se deduz do anterior quando h +» (27). 


3.1.0.2.3—Considera-se exemplificado o método, 
mas por ser de interesse nas aplicações, apresen- 
taremos sem dedução (28) devida a limitação de 
espaço as soluções que se obteriam no caso do 
cilindro e da esfera, de raio R homogéneos, 
isótropos e de propriedades constantes, a tempe- 
ratura inicial constante T1, h uniforme e o meio 
a T,, também constante. 


3.1.0,2.3.1 — No caso do cilindro, nas condições 
anteriores obtém-se para a equação de Helmoltz 
e em coordenadas cilíndricas. 
dg 


d'o 1 á 
E e E e) 
de r dr É 


3.30 


que é a equação de Bessel da ordem zero. 
Os valores próprios qk são as raízes (24) da 
equação transcendente. 
3.31 (gR)JN(gR)=HJo(gR) 
Joe Ji são as funções de Bessel de primeira 
espécie e de ordem zero e ordem um respectiva- 
mente. 


Heisler [16], etc. mas apenas as últimas abrangem a 
gama de valores que se encontram nas aplicações; as 
outras abrangem com certa precisão a curta região que 
englobam. Os de maisinteresse encontram-se na AEIST 
— Secção fotográfica para uso dos alunos de Caldeiras 
e Permutadores de calor. 

(27) É cómodo neste caso o tabelamento de 329 
feito para x==o por Olson e Schultz |16), 

(28) Pode consultar-se qualquer das obras já citadas 


[E, 2, 4, 18]. 
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A solução é dada por: 


3,32 DE a 
TT — To 
RR A 
SS 2 J(t) a E  ( =] 
== +» | Pk R 


Sendo pru=qk*R 

A série é também rápidamente convergente 
mesmo para valores relativamente pequenos do 
número de Fourier. 


3.1.0.2.3.2 — No caso da esfera, obter-se-ia: 


d? » 2 dp» 


3.33 
dr? r dr 


E os valores próprios q« são as raízes (24) da 
equação transcendente 


3.34 (1—-HR) = (gR) cotg (qR) 
e a solução é dada por 
gm TD. 
Ti— To mc 
- u um (o E) 
= E Sen UV COS | o E 
k=1 Vs — Sen v COS “k ve — 


sendo “wu = qr R 


3.1.0.3 — Tem interesse referir ainda a possi- 
bilidade de simplificação que por vezes existe na 
solução dos problemas quando as condições no 
contorno assumem certas formas particulares que 
permitem desdobrar o problema geral em pro- 
blemas mais simples. Tal é o caso, por exemplo, 
de problemas a três dimensões, em que a solu- 
ção se pode obter pelo produto de soluções de 
três problemas a uma dimensão, e que pelo seu 
interesse prático a seguir se exemplifica. 


3.1.0.3.1 — Suponha-se o caso da condução de 
calor num paralelipipedo rectângulo, em que o 
problema se exprime por: 

Determinar T(P,t) tal que: 


ii EE qo Tão A 
d xi dx) d Xs a dt 
Parat>0ea<xi<hbi 
a<x< b; 
a << x3< bs 


sujeito a 


dT 

3.37 — =H, 1 parax=aet>0 
d xi 
(1=1,2,3] 

ea 

) 

3.35 — =-—H; parax=b; t>O 
x 


sendo H; e H;! constantes (> 0) 

e tal que no instante inicial 

3.39 T(x1,x2, x3,0) =T; (x1,0) Ta (x2,0) Ts (x3,0) 
t>D0 


Neste caso o problema pode desdobrar-se nos 
três problemas seguintes a uma dimensão: 
Determinar T; (x;, t) 


tal que 

9º T; ) T; 

340 —i="Áparat>Oealx<h 
d Xi dt 
) T 

341 —— =H Tt>D0ex=a 
d xi 
dT 

3.42 i=-HTtDOex=b; 
UV Xi 


sendo no instante inicial 


Ti (x,0)=T(x,0) t=0 a<lx<b; 


com efeito, se Ti (xi t), Ta (x2t) e T; (xs t) são 
as soluções de cada um dos problemas a uma 
dimensão anteriores (i== 1,2, 3), a solução do 
problema a três dimensões (3.36 sujeito a 3.37, 
3.38 e 3.39) é 


3.43 T=T(xit).To(x2t) Ts (xs,t) 
como fácilmente se verifica substituindo em 
3.36, 3.37, 3.38 e 3.39. 


O interesse prático da simplificação é sobre- 
tudo manifesto quando a temperatura inicial é 
uniforme no paralelipípedo, porque então é ime- 
diato escrever 3.39 (29), Este é o caso que muitas 
vezes surge nas aplicações, e que aumenta a uti- 
lidade da representação gráfica adimensional das 
soluções a uma dimensão já anteriormente citada. 

O método estende-se imediatamente ao caso 
do cilindro finito, e a muitos outros/1). 

3.1.0.3.2 — Outros casos de simplificação re- 
sultam imediatamente do carácter linear da 
equação do calor, o qual permite muitas vezes 


(29) Note-se que se a temperatura inicial se exprime 
por uma série tripla de Fourier se obtém imediata- 
mente a forma 3.39. 


obter soluções somando soluções obtidas de pro- 
blemas mais simples. É o caso por exemplo de 
problemas com condições no contorno funções 
arbitrárias do tempo e outros. 

Como alguns exemplos se tornam evidentes 
usando o cálculo operacional citar-se-ão quando 
abordarmos futuramente este aspecto. Exempli- 
fica-se no entanto para um caso mais geral e 
complexo do que o tratado até agora, o qual 
permite estender o método de separação de va- 
riáveis a um modelo mais geral do que o consi- 
derado pois Q(P) 0 e as condições no con- 
torno não são homogêéneas. 


3.1.1 — Consideremos agora o caso geral de 
ser Q(P) + 0, isto é, do problema se exprimir 


por: 
Determinar T(P,t) 
Satisfazendo 
Prod CA qa — Qt) PemDet>o 
* dt eo 
sujeito a 


3.44 M (PS + HIP T=N(P) PemBet>O 


(com H (P) não idênticamente nulo) 


e T(P, 0) = F(P) Pem Det=0 


Neste caso, embora lineares em T tanto a 
equação como as condições nos limites não são 
homogéneas. 

O processo de resolução consiste na subdivi- 
são do problema nos dois mais simples seguintes. 

Determinar Z (P) tal que 


3.45 ta do ML sô P em D 
o 
M SÊ -HZ=eN P em B 
dn 
e W (P, t) tal que 
v? AE ug PemDet>o 
a dt 
M já +HW= Pem B et >0 
dn 


W=F(P)-Z(P) PemDet=o 


como facilmente se verifica, a solução do pro- 


blema proposto é dada por 
3.46 T(PG)=Z(P+AW(P,t) 
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O problema foi deste modo reduzido a um 
problema de regime estacionário — determinação 
de Z (P) —e a um de regime variável — W (P, t). 
Obtido Z (P), a determinação de W (P,t) corres- 
ponde ao problema já tratado em 3.70, que se 
sabe resolver. 

A determinação de Z que é um problema re- 
lativamente mais simples tratar-se-á posterior- 
mente. Aliás, como a uma dimensão é imediato, 
e este é o caso frequentemente necessário nas 
aplicações, as soluções anteriormente obtidas 
podem estender-se imediatamente a este caso. 


3.1.1.0 — Se H=0, 3.44 que corresponde a 
impor um fluxo térmico de valor determinado 
no contorno, o problema pode não ter solução 
no âmbito do método anterior [17]. Neste caso 
a solução pode obter-se somando a 3.46 um 
termo À t, sendo À uma constante dada por 


Das | fjnwas+ f om av] 
pCcv B D 


que na determinação de Z se deve somar ao se- 
gundo membro de 3.45. 


3.1.2 — Conclusão 


Dão-se seguidamente algumas indicações quanto 
à utilidade do método nos casos de aplicação. 
Algumas restrições foram explicitamente indi- 
cadas na formulação do modelo. Outras refe- 
rem-se a objecções de carácter prático. 

Assim, quanto ao sólido, este deve: 


a) Ser homogéneo (e de propriedades tér- 
micas constantes). 


b) Ser limitado (B) por superfícies coor- 
denadas. 


c) Ser limitado relativamente às coordena- 
das segundo as quais há variação de 
temperatura. 


Quanto às soluções obtidas: 


d) São úteis quando o número de Fourier 
(t =>) não é muito pequeno — caso 
contrário as séries convergem lentamente. 


e) Quando o método é aplicável são quase 
sempre as mais fáceis de obter. 
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Quanto a auxiliares de cálculo: 


As tabelas de raízes, ou tabelamento gráfico 
das soluções foram indicados já — estas existem 
para uma dimensão e usam-se para duas ou três 
dimensões quando o produto de soluções é pos- 
sível. As tabelas de raízes existem para as for- 
mas simples — paralelipípedo rectângulo, cilin- 
dro, esfera e cone [1). 

A bibliografia mais importante é a que se 
indica, a qual se ocupa exclusivamente, das for- 
mas geométricas simples — placa, paralelipípedo 
rectângulo, cilindro e esfera, com as restrições 
já indicadas, e quase sempre considerando a tem- 
peratura inicial uniforme [2, 3, 4, 13,14, 17, 18]; 
é excepção [1| que considera ainda o sólido 
limitado por superfícies cónicas, e casos de tem- 
peratura inicial não uniforme — é uma obra clás- 
sica em matemáticas aplicadas. 
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RESUMO 


Numa sucinta introdução procura-se definir correctamente o problema da condução do 
calor em sólidos invocando a Termodinâmica dos Processos Irreversíveis. Expõe-se seguida- 
mente o método de separação de variáveis de que se faz um resumo crítico relativo às van- 
tagens e limitações na sua utilização. 


SUMMARY 


The heat conduction in solids is treated in a brief introduction, in the scope of Irre- 
versible Termodinamics. Separation of variables method is then exposed, and its advanta- 
gens and limitations discussed. 
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Curso Adiantado de Engenharia Estrutural 


Faculdade de Engenharia — Universidade de Southampton 


Convidam-se todos os estudantes, devidamente qualificados, a inscreverem-se num curso de 
um ano de Engenharia Estrutural. Os candidatos, que terão de ser formados em Engenharia ou 
Matemática, ou possuirem sólidos conhecimentos destes ramos da Ciência em que desejem inscre- 
ver-se, deverão dirigir-se, o mais rápidamente possível, à Secretaria da Faculdade de Engenharia 
da Universidade de Southampton. 


O curso abrangerá os seguintes assuntos: 


Análise de estruturas elásticas lineares 

Teoria e prática do cálculo à rotura aplicada a estruturas de aço e de betão 
Estabilidade e vibração de estruturas 

Teoria das placas e das cascas 

Análise experimental de tensões 

Cálculo de fundações 

Teoria da elasticidade 

Métodos variacionais 

Uso de computadores electrónicos. 


assim como conferências sobre a Matemática relacionada com a teoria das estruturas. 


Fazem também parte do curso o estudo de investigação experimental e cálculo de 
projectos. 


Os estudantes já formados serão submetidos a um exame em Junho e devem apresentar uma tese 
no fim de Setembro, sendo-lhes conferido o grau de M. Sc. (Eng.) em caso de aprovação. 

Aos que frequentarem o curso com aproveitamento será atribuído um Diploma da Uni- 
versidade. 

As inscrições para o curso e os exames serão de 75 libras, excluindo alojamento. 

O Department of Scientific and Industrial Research reconhece o curso e oferece bolsas 
de estudo. 
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Moderno aparelho equi- 
pado com um sistema 
prismático pendular, co- 
locando automáticamente 
a linha de mira na sua 
posição horizontal. 

A simplicidade e robustez 
do processo utilizado e 
o seu grande rendimento 
de trabalho, tornam este 


aparelho excepcional 
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CARACTERISTICAS 


Luneta — Com focagem interna, óptica tratada, imagem direita, Comprimento 200 mm, 
Ampliação 20 x. abertura da objectiva 25 mm. Distância mínima de mira 2,5 m. 


Constante de adição zero, Constante ds multiplicação 100 


Limbo horizontal — Diâmetro 100 mm. Graduação 1º — | g 
Leitura directa 1º— 1 g Leitura por avaliação 6” — 10º 


Sensibilidade do nivel esférico: 8 
Tripé extensível, caixa metálica com pequenos acessórios 
Pesos: aparelho 2 kg, caixa 2,400 kg; tripé 5,200 kg 


REPRESENTANTES 


PAPELARIA FERNANDES — S.A. R.L. 


L. DO RATO, 13 — TEL. 68213139 — R. DO OURO, 145 — TEL 328361 / 3260 66 
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SUNWAY = Os Estores de concepção ultramoderna 


que graças às suas reconhecidas características de segurança, 
eficiência e poder decorativo são recomendados para 
qualquer ambiente. 


Sistema revolucionário de ligação das lâminas 

por meio de fios de Terylene extremamente finos 

e de extraordinária resistência, alta qualidade de todas 
as partes componentes, funcionamento silencioso 


e facilidade de limpeza — são as características principais 
dos Estores «Sunway». 


ESTORES SUNWAY 


(DIVISÃO DA BELL'S E CA, LOA, 


AVENIDA GENERAL ROÇADAS, 68-C « 68-00 
*qLes, B3 54 35 Tereo, SUNWAaN 
LISBOA 1 


ii 


C.D.U. 624.138 


NOTAS SOBRE A ESTABILIZAÇÃO DE SOLOS* 


I — GENERALIDADES SOBRE A ESTABILI- 
ZAÇÃO DE TERRENOS 


1 — ECONOMIA 


a) — Custo dos pavimentos 


O custe total da construção dum quilómetro 
de estrada pode considerar-se dividido nas seguin- 
tes parcelas: —terraplenagens e obras de arte, 
pavimento e obras acessórias. 

Em Portugal estes custos variam dentro dos 
limites apresentados no quadro 1. 


QUADRO I 


Custo de um km de estrada 
(Milhares de escudos) 


Terraplenagens Obras 
e obras de arte Pavimento acessorias Total 
Estradas mu- | 
nicipais | 
(Largura: 5m) 80 150 50 280 
Estradas na- 
cionais 
(Largura:7 m) 300 350 200 850 
Auto-estradas 
(Largura : 
2X7m=14m) | 2500 1500 2500 | 6500 
Custo relativo do pavimento 
Estradas municipais ... 54% 
Estradas nacionais .... 41% 
Auto-estradas . ...... 23 “Jo 


Quer dizer, o custo do pavimento é tanto mais 
importante quanto menos importante é a estrada. 
Daí o grande interesse do conveniente estudo 
dos pavimentos nas estradas de baixo custo. 

O custo da conservação, que é de muito difícil 
avaliação, é dum modo geral tanto maior quanto 
mais barato é o pavimento. Nas auto-estradas 
esse custo será relativamente pequeno, mas nas 


por ÚLPIO NASCIMENTO 


Chefe do Serviço de Geotecnia do Labo- 
ratório Nacional de Engenharia Civil 


estradas de baixo custo a conservação poderá 
atingir 14 ou 1/3 do custo da construção. 

No quadro II apresentam-se as ordens de gran- 
deza dos custos em Portugal metropolitano de 


QUADRO II 


Custo em escudos por m? de diferentes camadas utiliza- 
das como fundação e camada de desgaste de pavimentos 


1 — Terra estabilizada com 15 cm de 
espessura depois do recalque (se- 
leccionada e homogeneizada). .. 

2 — Macadame com 15 cm de espes- 
sura depois do recalque (0,18 mº 
de brita a 45$00/mº) ........ 

3 — Solo-cimento com 15 cm de es- 
pessura depois do recalque (22 
kg de cimento a $70/kg). .... 

4 — Revestimento betuminoso sobre 
macadame (2,6 kg de betume a 
RERROUMES cavaco é ums é Reu à 

5 — Revestimento betuminoso sobre 
base de terra estabilizada 
(=a4 +11 de MCO a 2$00/1) 

6 — Tapete betuminoso com 4 cm de 
espessura (5,5 "/» de betume a 
1$70 kg + 0,6 kg de emulsão a 
1$70 kg em rega de colagem'.. 

7 — Idem com 7 cm de espessura 
(5,5 */o de betume + 0,6 kg de 
emulsão em rega de colagem)... 

8 — Calçada de paralelipípedos de 
granito (70 cubos a $25 cada + 
+ 0,12 mº de areia a 30$00/m?) 

9 — Laje de betão de cimento de 
15 cm de espessura (betão de 
350 kg de cimento a $70 kg)... 

10 — Laje de betão armado de 15 cm 
de espessura (3,8 kg de arma- 
dura a SSOO/kE) «isesucscs 


10$00 


15$00 


32$00 


13$00 


16$00 


25$00 


41$00 


36$00 


* Conferências realizadas em Madrid, a 9 e 10 de Janeiro de 1960, num curso de construção de estradas organi- 


zado pela Asociación Espafiola de la Carretera. 
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algumas camadas utilizadas como fundação é 
camada de desgaste. 


b) — Capacidade de tráfego dum pavimento 


Um pavimento submetido a um dado tráfego 
tem uma certa duração além da qual a sua ruína 
se acentua progressivamente manifestando-se 


por fendilhamento, ondulações, rodeiras, covas, crema- 
malheiras, etc. Tem portanto uma certa capacidade 


de tráfego que se poderá medir pelo número má- 
ximo de veículos a que pode dar passagem em 
boas condições. 

Sendo I a intensidade diária média de veículos 
durante um ano e D a duração útil da estrada, 
em anos, a capacidade do pavimento será dada 
pelo número total N de veículos que se obtém 
pela expressão 


N = 365 ba Th (1) 


Se designarmos por I* o valor médio de I ao 
longo de toda a duração D da estrada, à eq. (1) 
poderá dar-se a forma 


N=365 ID (2) 


Evidentemente que esta capacidade depende 
fortemente da natureza do tráfego que circula na 
estrada, principalmente da sua carga por roda, 
A experiência portuguesa permite no entanto 


fixar no quadro III algumas ordens de gran- 
deza de D, I e N. 


QUADRO III 
Capacidade de tráfego (N) de alguns pavimentos 
N=3651'D 
Tipos | xa ' y | 
P (anos) (velculos[dia)| (veiculos) 
Serpa Essas ne) id | 
1 — Terra estabilizada . | 5a 10 10 '2a4x<10 
2 — Macadame sem re- 
vestimento . .«. . .| 5a1o 100 2a4>(10º 
3 — Revestimento su- 
perficial sobre terra 
estabilizada, . . .|! 5a1lo 1000 2a4 x 10º 
4 — Revestimento su- 
perficial sobre ma- 
cadame +... .| Saio 1500 3a 6X 10º 
5 — Tapete betuminoso |10a 20 3000 1a2><10 
6 — Laje de betão de 
cimento . « . .« .« |15a30 3000 


PR 
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c) — Custo efectivo dos pavimentos 


O valor económico dum pavimento é medido 
pelo serviço que ele presta. Uma maneira gené- 
rico de medir esse serviço será pelo número de 
veículos que passam sobre um quilómetro da 
respectiva estrada, isto é, pelo número de veículos. 
km que ele pode produzir. O valor económico 
dum pavimento será pois o do seu custo num 
quilómetro de estrada, dado pela expressão 


P=10007p1 (3) 


sendo 


7 — um coeficiente, para dar conta do custo da 

conservação, variando entre 1 e 1,3; 
p—custo de um m” do pavimento (quadro II) 
l—a largura da estrada em metros, 


dividido pela sua capacidade, expressa pelo nú- 
mero máximo de veículos que por ele pode passar. 

Obter-se-á assim o custo/veículo.km (c) do 
pavimento pela expressão 


P 1000 | 
a E EPE RA (4) 
N N 
QUADRO IV 
Custo efectivo por veículo. km de alguns pavimentos 
po SER (y=11|=7,5m) 


Custo | Capacidade Custo efectivo 


Tipo caca diê N (escudos por 
DOE ni?) (veículos) veículo, km) 


zada .« . . .« .|10a20 |2a4> 10" 2a 8 
2 — Macadame sem 
revestimento. . 
3 -—Revestimento 
superficial sobre 
terra estabili- 
gada ww wo 
4 — Revestimento 
superficial sobre 
macadame. . 
5 — Tapete betumi- 
| Dodo: » » & x 
6 — Laje de betão de 
cimento , . .« 


15a30 /22a45<10º| 03a1, 


25a35 22a4X10*|, 0,054 0,13 


30a 45 326 10º| 0,042 0,11 


40 a 70 


1— Terra estabili- 
, a2>< 10º |0,015a 0,09 


80 a 100 2 a3> 10' | 0,02 a 0,04 


No quadro IV apresentam-se algumas ordens 
de grandeza deste custo efectivo deduzidas dos 
quadros II e III, supondo uma estrada com 


|==7,5 m de largura e não entrando em conside- 
ração com os encargos de conservação (y =1). Por 
aí se verifica que os pavimentos mais caros são 
também os mais económicos, isto é, aqueles cujo 
custo efectivo por veículo. km é mais baixo; e 
verifica-se também que os pavimentos mais ba- 
ratos são os menos económicos. As diferenças 
podem ir desde 2 centavos/veículo . km, nos me- 
lhores pavimentos betuminosos ou de cimento, 
até 8 escudos/veículo.k numa estrada de terra; 
quer dizer a relação entre estes custos extremos 
pode ir até 400 (8/0,02). 


d) — Planeamento por fases 


Poderia então parecer que as estradas de baixo 
custo não têm interesse económico. Há porém 
que entrar também em consideração com o as- 
pecto financeiro do problema. Com efeito se o 
pavimento tem elevada capacidade exige um 
grande investimento inicial e para isso será ne- 
cessário dispor de dinheiro que só poderá ser 
amortizado ao fim dum grande número de anos. 
Haverá então que somar ao custo do pavimento, 
os encargos com os juros do financiamento; e 
nestas condições já o confronto com os pavi- 
mentos de baixo custo pode ser favorável a 


estes. 
Mas há ainda outras considerações de muita 


los) 


tcu 


* 


apacidade de trafego N (ve 


Ç 


10 20 30 40 
Duração D (anos) 


Fig. 1 — Esquema da evolução do tráfego em duas 
estradas, uma de pouca (1) e outra de muito tráfego (2), 
e da construção por fases no caso (1). 


importância que dizem respeito à previsão da 
evolução do tráfego. 


O serviço que um quilómetro de estrada vai pres- 
tando ao longo do tempo é medido pelo número 
de veículos.km dado pela eg. (1) o qual pode 
evoluir, por exemplo, conforme se representa na 
fig. 1, segundo as curvas (1), para estradas de 
pequeno tráfego, ou (2) para estradas de tráfego 
mais intenso. 

Um bom pavimento, cujo custo efectivo é, 
como vimos, o mais baixo, levaria trinta ou qua- 
renta anos a gastar toda a sua capacidade de 
tráfego N, se este crescer muito lentamente con- 
forme se exemplifica com a curva (1); enquanto 
que se gastaria apenas em dez ou vinte anos, se 
o tráfego crescer rápidamente conforme a curva 
(2). No caso (1) o financiamento seria a muito 
longo prazo o que terá os seguintes inconve- 
nientes : 


1.º — Grandes encargos com juros, o que vai 
J 
elevar o custo efectivo do pavimento. 


2.º — Exigir uma grande confiança na previsão 
da evolução do tráfego ao longo de várias 
dezenas de anos, o que constitui proble- 
ma que ainda está longe de solução 
satisfatória. 


Por estas razões é preferível fazer uma cons- 
trução por fases de capacidades Nº e N” (fig. 1) e 
durações D' e D” de modo que a estrada se vá 
adaptando por degraus ao crescimento do tráfego. 
O custo efectivo desta solução poderá ser maior; 
no entanto tem as seguintes vantagens : 


1.4 — Reduz a duração de cada fase a um prazo 
dentro do qual são razoáveis as previsões 
do tráfego e os investimentos financeiros 
respectivos. 


2.2 — Permite que a estrada, pela contribuição 
de cada uma das suas fases para a econo- 
mia do país, vá ela própria ganhando o 
dinheiro que permitirá construir a fase 
seguinte. 


É nesta construção por fases, principalmente 
nos territórios do ultramar em desenvolvimento, 
que as pavimentações com terras estabilizadas 
têm mais interesse. 

Todavia para se poder fazer o planeamento é 
necessário conhecer em relação a cada fase a 
respectiva capacidade N e a duração D. Sabemos 
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nós calcular tais valores com rigor suficiente? 
A resposta é infelizmente negativa. 

Com efeito o cálculo da capacidade de tráfego N 
para um dado pavimento constitui um impor- 
tante problema de dimensionamento que está 
longe de estar resolvido, sobretudo para pavi. 
mentos que saiam fora das condições normais, 
como é geralmente o caso da construção por fases. 
Por outro lado o cálculo da duração D constitui 
um problema de previsão, baseado nas estatis- 
ticas de tráfego ou em planos de desenvolvimento 
económico, que também não se pode considerar 
resolvido na generalidade dos casos. 

Mas então, se não conhecemos com suficiente 
rigor os parâmetros em que se baseia o planea- 
mento qual é o seu interesse? O seu interesse é 
o de pôr o problema em equação em bases quan- 
titativas, e isso é o primeiro passo para a sua 
solução. 


2 — ESTUDOS, ENSAIOS E PROJECTO 


O projecto duma estabilização deve ter por 
base os seguintes estudos: 


— Reconhecimento geológico dos terrenos ao 
longo da estrada e até uma distância desta 
dentro da qual ainda se justifique a utiliza- 
ção económica dos respetivos materiais. 

— Prospecção, amostragem, identificação e clas- 
sificação dos terrenos com interesse para a 
obra. 

— Ensaio dos materiais tendo como base a sua 
utilização nas possíveis soluções de pavi- 
mentação consideradas em anteprojecto como 
mais económicas do ponto de vista, quer da 
melhor utilização dos materiais da região, 
quer dos meios de trabalho disponíveis para 
a construção. 

— Dimensionamento. Definida a constituição 
da camada de terra estabilizada e caracteri- 
zada a respectiva resistência, pelo C. B. R. 
por exemplo, haverá que dimensioná-la, isto 
é, calcular a sua espessura. 


Vejamos cada um destes pontos com mais 
pormenor: 


a) — Reconhecimento geológico 
O reconhecimento geológico é recomendável 
porque define e localiza as principais formações 
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da região réferindo-as à geologia geral do país 
o que apresenta as seguintes vantagens: 


1.2 — Permite tirar partido de elementos infor- 
mativos que já existam sobre outras ocorrências 
da mesma formação. E, reciprocamente, facilita a 
futura utilização noutros trabalhos das informa- 
ções que se reunam no trabalho em questão. 

2.2 — Permite orientar com maior eficiência e 
economia a amostragem de pormenor para os 
estudos geotécnicos. 

O reconhecimento geológico deve ser prece- 
dido da reunião da bibliografia sobre a geologia 
da região (cartas geológicas, monografias, etc.) e 
de outros elementos informativos disponíveis, 
tais como relatórios de sondagens feitas na re- 
gião mesmo para outros fins. Com base nestes 
elementos e nas cartas topográficas existentes 
faz-se um reconhecimento no terreno, localizando 
e identificando, se possível, as formações dos 
afloramentos das trincheiras, dos barrancos, etc. 
Deste modo é possível estabelecer uma primeira 
planta onde se delimitam as zonas em que afloram 
as formações ou os complexos geológicos identifi- 
cados. Pode mesmo ser possível esboçar alguns 
cortes geológicos de muito interesse para o es- 
tudo ulterior. 

Com base no relatório do reconhecimento 
geológico estabelecer-se-á o plano das sondagens, 
dos poços e das valas de prospecção que consti- 
tuem os trabalhos da fase seguinte. 


b) — Prospecção, amostragem, identificação e classifi- 
cação dos terrenos 


Os principais objectivos destes trabalhos são: 


1.º — Confirmar e completar as informações 
reunidas no reconnecimento geológico, 
delimitando em superfície e em profun- 
didade as principais formações geológi- 
cas com interesse para obra. 

2.º — Fazer a descrição da composição litoló- 
gica e da textura e estrutura dos mate- 
riais com base no exame visual seguin- 
do-se por exemplo a orientação da 
norma AASHO — T 86-49. 

3.º — Fazer a amostragem, identificação e clas- 
sificação geotécnica dos materiais com 
base na granulometria e limites de Atter- 
berg seguindo-se por exemplo a especi- 
ficação AASHO M 145-49, 


Esta amostragem e identificação devem 
ser orientadas de maneira a poder fazer-se 
uma ideia da heterogeneidade dos mate- 
riais e das respectivas propriedades 
médias. 


O poço ou a vala que permitem o exame di- 
recto dos materiais «in situ» e a sua amostragem 
nas melhores condições são os métodos mais 
recomendáveis de prospecção. Em certas circuns- 
tâncias (ocorrência de água, grandes profun- 
didades, etc.) terá que se recorrer aos furos de 
sondagem. 

Uma questão fundamental é a do afastamento 
entre si e da profundidade dos poços e sonda- 
gens de prospecção. Não me demorarei neste 
problema porque alongaria demasiadamente este 
trabalho. 


c) — Ensaios 


Os ensaios que correntemente se realizam so- 
bre os materiais têm como principal objectivo 
uma primeira solução em face dás especificações 
adoptadas. Quer dizer uma primeira avaliação do 
seu valor em face duma experiência prática que 
é geralmente divulgada e oficializada por meio 
de normas e especificações. Entre nós, como em 
muitos outros países, tomam-se como ponto de 
partida as especificações americanas que se vão 
adaptando e oficializando à medida que as respec- 
tivas técnicas se vão desenvolvendo. 

Há porém que ter presente que tais especifi- 
cações foram elaboradas para atender a circuns- 
tâncias, geológicas, climáticas e económicas que 
podem ser muito diferentes das do nosso país. 
Portanto, mesmo nos casos em que elas sejam 
aplicáveis, há que manter uma certa reserva 
enquanto não houver suficiente experiência 
acumulada através das respectivas aplicações 
Quando tal não sucede a melhor orientação será 
a de adquirir essa experiência o mais rápida e 
econômicamente possível. Para isso impõe-se uma 
intima colaboração entre as organizações de in- 
vestigação e as organizações de construção e 
conservação de estradas, no sentido de se cons- 
truirem numerosos e bem distribuídos troços ex- 
perimentais, aproveitando para o efeito os pró- 
prios programas de construção de rede nacional, 
nos quais se comprovem em condições reais os 
projectos de especificações definidas pela investi- 
gação laboratorial. 


Adiante, a propósito de cada tipo de estabili- 
zação, referirei com pormenor as respectivas es- 
pecificações e ensaios. 


d) — Dimensionamento 


As duas modalidades de planeamento repre- 
sentadas pelas curvas (1) e (2) da fig. 1 dão 
também os dois casos mais correntes de utiliza- 
ção de camadas de terra estabilizada: 


1.º — A camada está integrada num pavimento 
cuja duração se supõe normal de acordo 
com os critérios correntes de dimensio- 
namento (C. B. R.). 

2.º — A camada ou o pavimento a que ela per- 
tence constituem uma fase cuja duração 
só será normal sob tráfego muito redu- 
zido; ou então, se o tráfego for normal, 
terá uma duração muito breve. 


Com efeito, no que diz respeito ao tráfego, 
os parâmetros fundamentais do dimensionamento 
são a carga da roda representativa do tráfego e 
O número de repetições dessa carga. Em primeira 
aproximação esse número de repetições pode 
considerar-se proporcional ao número total (N) 
de veículos que passam na estrada, isto é, à sua 
capacidade de tráfego dada pela eg. (2). 

Os métodos de dimensionamento corrente- 
mente adoptados, designadamente o método 
C. B. R., não tomam em linha de conta duma 
forma explícita esta capacidade de tráfego, pois 
apenas se referem à carga da roda. Subentende- 
-se porém que tais métodos se destinam a es- 
tradas de tráfego normal e com uma duração 
também normal. 

Apenas para referir ordens de grandeza con- 
sideraremos normal uma duração de 15 a 20 anos 
com uma intensidade de 1000 a 2000 veículos 
por dia; teremos assim que a capacidade de tra- 
fego das estradas mencionadas pelo C. B. R. é 
de 10º a 10” veículos. 


Ncer = 10º a 10” veículos 


Consideremos o exemplo apresentado na fig. 2 
de um pavimento dimensionado pelo C. B. R. 
para rodas de 5,5 ton (12 000 libras) sobre um 
terreno de CBR = 5. Tem que ter 45 cm de 
espessura total, a qual pode ser realizada por um 
tapete betuminoso de 7 cm, uma base de maca- 
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me de 13 cm e uma sub-base de terra estabili- 
zada com CBR = 20, OC. B.R. da base tem 
que ser elevado, pelo menos 80, para que, de 
acordo com o ábaco, a espessura do tapete possa 


Mg +“ 
o e 
E Bend. data 
ER 


45 cm 


25 


CAMADA DE DESGASTE 


BASE .. MACADAME 


” - , e PRESO E 3% 
DMI CBR=20." SUB-BASE - TERRA ESTABILIZADA 


TERRENO DE FUNDAÇÃO 


a sua espessura começa a ser diminuida à custa 
da acumulação de materiais nas bermas. 

Se dermos à camada apenas os 25 cm que lhe 
virão a corresponder de futuro ela será desagre- 


= 
So 


ESPESSURA (cm) 


Fig. 2 — Exemplo de dimensionamento de pavimentos pelo C. B.R. 


ser apenas de 7 cm. Um pavimento nestas con- 
dições terá portanto uma capacidade de 10º a 107 
veículos e durará 15 a 20 anos submetido a um 
tráfego médio de 1000 a 2000 veículos por dia. 

No entanto se a estrada nos primeiros anos 
tiver um tráfego muito pequeno, poderá justifi- 
car-se a construção por fases como exemplificá- 
mos na fig. 1 (curva 2). 

Suponhamos que são as três seguintes fases: 


1.º fase — Sub-base de terra estabilizada com 
revestimento superficial betuminoso. 

2.º fase — Regularização e recarga com maca- 
dame e revestimento betuminoso. 

3.º fase — Regularização e recarga com tapete 
betuminoso. 


Voltemos porém à fig. 2 e consideremos o pro- 
blema concreto do dimensionamento duma ca- 
mada de terra estabilizada quando constitui a 
1.º fase. 

Se dermos à camada toda a espessura de 45 cm 
que o terreno de fundação requer não se satisfaz 
o método do C. B. R. nos 20 cm de superfície 
visto que aí são requeridos materiais com C.B.R. 
superior a 20 e tanto maiores quanto mais pró- 
ximo da superfície eles se encontram (fig. 2). 
A terra estabilizada será portanto rapidamente 
desagregada dando lugar a rodeiras, lamas, poei- 
ras, cremalheiras, etc. E como consequência disso 
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gada na mesma, mas com a grande agravante de 
tal espessura ser insuficiente para o €. B. R. do 
terreno de fundação, o que implica uma redução 
apreciável da sua capacidade de tráfego e por- 
tanto a eventual destruição de todo o pavimento 
desde a sua base, impedindo pois o seu aprovei- 
tamento para as fases futuras. 

A solução a adoptar em cada caso, que depende 
evidentemente das respectivas circunstâncias, 
pode localizar-se entre as duas seguintes situa- 
ções extremas: 


1.º — A camada tem apenas a espessura que 
lhe corresponderá como sub-base do futuro pavi- 
mento depois de construídas todas as fases 
(25 cm no exemplo da fig. 2). 

2.º — A camada tem, pelo menos, a espessura 
requerida pelo C. B. R. da fundação (45 cm no 
exemplo da fig. 2). 

A construção das fases seguintes conduzirá 
portanto a um pavimento sobredimensionado. 

Todavia, em qualquer dos casos, como não 
conhecemos nem a capacidade N nem o número 
de anos D que ela levará a ser atingida, não 
sabemos também «à priori» qual será o momento 
em que será necessária a construção da 2.º fase. 
Aliás essa necessidade da 2.º fase pode ser con- 
siderávelmente retardada se houver uma cuida- 
dosa conservação e mesmo trabalhos de pro- 
tecção. 


A boa orientação parece ser a seguinte: 


1.º — Manter permanentemente a estrada sob 
cuidadosa e sistemática observação. 

2.º — Realizar sem demora os trabalhos de 
conservação e de protecção sempre que 
necessários. Uma despesa feita a tempo 
pode poupar futuras despesas muito mais 
avultadas. 

3.º — Construir a fase seguinte antes que a 
ruína seja tão profunda que a torne inu- 
tilizável como fundação, sendo portanto 
necessário levantá-la e refazer todo o 
pavimento desde o terreno de fundação. 


II — ESTABILIZAÇÃO MECÂNICA 
1 — Especificações 


A estabilização mecânica é aquela que se obtém 
apenas por operações de carácter mecânico, tais 
como lavagem, crivagem, britagem, mistura, hamogenei- 
zação e compactação dos materiais. 

Pretende-se melhorar a sua resistência às solici- 
tações do tráfego e diminuir a susceptibilidade 
dessa resistência às variações de humidade. 

A resistência mecânica depende principalmente 
da compacidade : - quanto mais elevado for o 
peso específico aparente dum terreno maior será 
a sua resistência. E essa compacidade depende 
principalmente da granulometria. 

A susceptibilidade à humidade é devida à ex- 
pansibilidade dos finos do terreno provocada pela 
absorção de água, pois esse aumento de volume, 
diminuindo a capacidade, diminui a resistência. 
Analogamente a excessiva retracção dos finos por 
secagem pode também provocar grande queda 
da resistência e a consequente desagregação do 
material. No primeiro caso o terreno forma lama; 
no segundo caso forma poeiras. A estabilização 
tem que evitar uma coisa e outra. Vejamos como. 


a) — Granulometria 


Quanto à granulometria, está averiguado serem 
aquelas que se aproximam das do tipo Fuller ou 
Talbot as que conduzem ao maior peso específico 
aparente. No entanto, para melhorar a coesão, 
a quantidade passada no peneiro n.º 200 tem 
que ser aumentada. Do ponto de vista das apli- 
cações práticas a especificação ASTM D 1241- 
-55 T define vários fusos granulométricos para 


dois tipos de materiais: o tipo 1 com 4 varieda- 
des (A, BC e D) de materiais mais grosseiros 
que poderemos designar por betão de solo e o 
tipo II, com duas variedades de materiais mais 
finos (E e F), que poderemos designar por arga- 
massa de solo, 

Na fig. 3 apresentam-se esses fusos para as 
quatro variedades de betões de solo (A,B, Ce D) 
e na fig. 4 as duas variedades de argamassas de 
solo (E e F). Além das limitações destes fusos, 
impõe ainda a mesma especificação ASTM que 
a quantidade de finos passando no peneiro 
n.º 200 (0,074 mm) não exceda 2/3 da quantidade 
que passa no peneiro n.º 40 (0,42 mm). No 
entanto, se se tratar duma camada de desgaste, 
essa quantidade de finos passados no peneiro 
n.º 200 não pode ser inferior a 8º/. 


b) — Plasticidade 


O comportamento perante a água é contro- 
lado por meio dos limites de Atterberg. Impõe- 
-se que a fracção passando no peneiro n.º 40 
tenha 

LL< 25 e IP< 6 


No entanto, tratando-se duma camada de des- 
gaste, deverá subir-se o limite superior e limitar 
inferiormente a plasticidade para evitar a desa- 
gregação e as poeiras 


LL; *<E.:35 [Ê==4 a 9 


c) — Dureza do agregado grosso 


Quanto à dureza do agregado exige-se ainda 
que o material retido no peneiro n.º 10 (2 mm) 
tenha um desgaste na máquina de Los Angeles 
não superior a 50º/. 


d) — Resistência e compactação 


O método mais divulgado para a avaliação da 
resistência dum terreno para efeitos de dimen- 
sionamento de pavimentos é o do C.B. R. Dado 
um terreno satisfazendo às especificações de gra- 
nulometria e plasticidade referidas anteriormente, 
há pois que avaliar qual deverá ser o C. B.R. a 
adoptar no projecto tendo bem em atenção as 
circustâncias especiais e o equipamento com que 
se vai realizar a compactação que terá de garan- 
tir o referido €, B. R. 
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SEDIMENTAÇÃO PENEIROS Nº8 
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Fig. 3 — Especificação ASTM : D 1241-55 T para betões de solo 
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Fig. 4 — Idem para argamassa de solo 
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AÇO TOR 


O AÇO DE ALTA RESISTÊNCIA 
PARA BETÃO ARMADO COM 
MAIOR VENDA NO MUNDO 


TENSÃO DE SEGURANÇA 2 400 Kg /cm 2 
SEM GANCHOS 


Especifique 
AÇO TOR 
com o símbolo 


D 


SOABAL - SOCIEDADE DE AÇOS PARA BETÃO ARMADO, LDA. 
RUA JOAQUIM BONIFÁCIO, 2, 1.º — TELEF. 405 66 — LISBOA 


TÉCNICA — XIX 


FEITO q 4 
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e 


ferro fundido 

ferro fundido com grafite esferoidal 
aço vasado 

aço inoxidável 


DN 10 mm até 600 mm PN 6 a 40 kg/cm 


conforme as normas DIN e ASA 


VÁLVULAS DE SEGURANÇA, DE REDUÇÃO, 
ELECTROMAGNÉTICAS E VÁLVULAS ESPECIAIS 


SOCIEDADE PORTUGUESA DE VÁLVULAS, LDA. 


RUA DA ACADEMIA DAS CIÊNCIAS, 5 — LISBOA 2 
TE LERONE: PPCA 31710—5 LINHAS 
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Os principais factores de que depende a resis- 
tência do terreno avaliada pelo C. B. R. são a 
dimensão máxima da sua granulometria, O seu grau 
de compactação e a sua humidade. 

Examinemo-los sumariamente: 


1.º — Dimensão máxima 


Em relação ao 1.º factor, enquanto que as espe- 
cificações granulométricas consentem, como vi- 
mos, estabilizações com terrenos até 2”, o ensaio 
C. B. R. não pode ser feito com partículas maio- 
res do que 3/4”. Todo o material retido no 
peneiro de 3/4” tem que ser substituído por 
igual peso de material passado no 3/4 e retido 
no n.º 4. 

O material ensaiado em laboratório pode, em 
certos terrenos mais grosseiros, ter cerca de 40º/ 
dos seus elementos mais grossos (superiores 
a 3/4”) substituídos por partículas mais finas 
(inferiores a 3/4”). 


Por esta razão no Laboratório Nacional de 
Engenharia Civil está em estudo a possibilidade 
de ensaiar materiais sem limite de dimensões 
utilizando-se para o efeito moldes de grandes 
dimensões com diâmetros de 30, 46 e 60 cm 
(fig. 5). 

O diâmetro do molde terá que ser cerca de 
3 vezes o diâmetro da placa de carga e este, por 
sua vez, cerca de 3 vezes a maior dimensão das 
partículas do terreno. 

Um terreno com partículas de 2” (5 cm) exigirá 
portanto um molde de (3x3) x 5 = 45 cm de 
diâmetro e uma placa de carga de pelo menos 
15 cm de diâmetro. 

O resultado do ensaio é dado pelo módulo de 
resistência (3) obtido pela expressão 

asd 


A ea (5) 


sendo 


R — módulo de resistência (kg/cm?) 


Fig. 5 — Moldes com diâmetros de 15 (C.B.R.), 30, 46 e 60 cm para ensaio de terrenos e agregados 
com dimensões máximas de 19, 35, 55 e 75 mm, respectivamente 


Esta limitação resulta das dimensões do molde 
C. B. R. e do diâmetro do respectivo pistão em 
relação ao tamanho das partículas. A resistência 
do material é avaliada por defeito e portanto o 
erro que se comete está do lado da segurança. 

Esta limitação não permite porém tirar todo o 
partido que se poderia tirar de materiais mais 
grosseiros e de maior resistência. 


q — tensão exercida sobre a superfície de 
terreno (kg/cm?) 

d — diâmetro da placa de carga (cm) 

3 — assentamento da placa para a tensão 9 (cm) 


A relação do módulo de resistência com o 
C.B.R. é sensivelmente a seguinte: 


R=(30 2:30) :5€ GC. BR. (6) 
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QUADRO V 


Tipos de ensaio de compactação segundo as normas ASTM 


| Compactação ligeira 
ASTM D 698-58 T 
Pilão: 5,5 libras; Queda: 12” 


Molde de 
Material passado 4.584“ de 
no peneiro alturá mem mm Nº de 
Tipo camadas 
E O ne 
| diâmetro: 4” | A 3 
N.º 4 (4,76mm) RR aanas 
diâmetro: 6” | B 3 
diâmetro: 4” | €C | 3 
3/4” (19,1 mm) E O 
diâmetro: 6” | D | 3 


| Compactação pesada 


ASTM D1557-58 T 
Pilão: 10 libras; Queda: 18” 


N.º de | Nº de | N.º de 
pancadas / cam. A] camadas | pancadas / cam. 
25 A| [2 
56 B. 5 | 56 
o 25 | e 5 E E 25 

56 “ID 5 | 56 


2.º — Grau de compactação 


O grau de compactação pode ser dado em rela- 
ção ao ensaio de compactação ligeira ou ao ensaio 
de compactação pesada, geralmente designado por 
ensaio AASHO modificado. 

Recentemente a ASTM normalizou esses ensaios 
permitindo que tanto a compactação ligeira como 
a pesada sejam feitas no molde de 4” ou de 6” e 
ainda, com material passado no peneiro n.' 4 ou, 
no peneiro 3/4”. Há portanto oito tipos de com- 
pactação que se resumem no quadro V. 

A qual destes tipos deve ser referido o grau 
de compactação? 

As compactações dos tipos A e B (N.º 4) estarão 
mais adaptadas para os terrenos finos do tipo IJ, 
enquanto que as compactações dos tipos Ce D 
(3 4”) responderão melhor nos terrenos mais 
grosseiros do tipo 1. 

Quanto ao tamanho do molde de compactação, 
parece haver algumas vantagens no molde maior, 
principalmente para os terrenos mais grosseiros, 
Uma dessas vantagens será uma mais estreita 
relação entre os ensaios de compactação e os 
ensaios €. B.R. 

O tipo de compactação a adoptar, ligeira ou 
pesada, é, em princípio, indiferente visto que o 
ensaio apenas constitui um padrão de referência, 
não limitando em nada a compactação a realizar 
no campo. Se a compactação escolhida for a ligeira 
pode exigir-se por exemplo um grau de compac- 
tação não inferior a 100 “/o; se for a pesada bas- 
tara, possivelmente, exigir 90 "/o. Note-se porém 
que não há uma equivalência geral entre os graus 
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referidos a um e outro tipo de compactação: 
cada terreno terá a sua equivalência que os res- 
pectivos ensaios poderão estabelecer. 


3.º — Humidade 


O €. B.R. dum terreno desce em regra muito 
rapidamente quando a sua humidade aumenta. 
Por isso o ensaio é geralmente feito depois da 
amostra ter estado mergulhada em água durante 
quatro dias. Pretende-se assim reproduzir no la- 
boratório as piores condições em que o terreno 
se pode encontrar na estrada. 

Em certos casos, porém, em que haja boas con- 
dições de protecção e drenagem do terreno, bas- 
tante experiência na região ou estudos especiais 
suficientemente desenvolvidos, poderão tomar-se 
valores do C. B. R. correspondentes ao terreno 
com menor embebição. 


2 — Projecto duma estabilização mecânica 


O projecto duma estabilização mecânica deve 
estabelecer todas as características a que deverá 
satisfazer a obra e, além disso, definir a maneira 
pela qual se controlarão essas características. As 
principais características a definir são a homege- 
neidade, a granulometria, a plasticidade, O peso específico 
aparente e a espessura da camada a estabilizar ; e as 
operações respectivas serão a homogeneização, a 
correcção granulométrica, a correcção da plasticidade, e a 
compactação. Analisemo-las sumariamente: 


a) — Homogeneização 


Uma das principais características que se exige 


dum pavimento é a sua regularidade, isto é, o per- 
feito desempeno da sua superfície, e que ela ao 
longo do tempo se deforme igualmente em toda 
a sua extensão sem provocar covas e ondulações. 
Grandes variações de resistência de ponto para 
ponto do terreno provocam diferentes assenta- 
mentos perante o tráfego e portanto a irregulari- 
dade do pavimento. 

Que haja locais de elevada resistência pouco 
interessa, se a pequena distância há outros de 
baixa resistência. É mesmo preferível uma resis- 
tência média baixa mas homogénea do que uma 
mais elevada mas heterogênea. 

Como definir a homogeneidade ? 

O assentamento provocado por uma roda num 
pavimento é práticamente o resultado da defor- 
mação dum volume de terreno cujas dimensões 
são em diâmetro e altura 2 a 3 vezes o diâmetro 
da área carregada. Como a estrutura do pavi- 
mento distribui as cargas na fundação este diá- 
metro da área carregada é maior do que o diá- 
metro da área de contacto do pneu que anda por 
20 a 30 cm. 

O volume do terreno deformado tem pois a 
dimensão de 0,5 m a 1,0 m. 

Numa primeira análise muito grosseira pode 
admitir-se que os assentamentos são inversa- 
mente proporcionais ao C.B.R. e portanto es- 
crever. 

C 


— (em) (7) 
CBR 

C—É uma constante cujo valor depende da 
pressão do pneu e do diâmetro da área de con- 
tacto. 

O seu valor em ordem de grandeza poderá 
deduzir-se das expressões (5) e (6) 

sd 


C=———— em 
20 a 30 


para 7=5 kg/cm? e d =30 cm 


Entre pontos do terreno afastados de dx haverá 
um assentamento diferencial ou irregularidade do 
cuja expressão se obtém derivando (7) 


C 


(C.B.R.)* 


d(C.B.R.) (8) 


Daqui se tira que as diferenças de C.B.R. entre 
pontos afastados de dx não poderão ser supe- 
riores a 


d(C.B.R.) < dá (9) 


(C.B.R.)º 
E a 
para que as irregularidades entre os mesmos 
pontos afastados de dx não sejam superiores a dº. 
Exemplifiquemos : 
Admitamos que dx ==1 m e que não queremos 
que dº seja maior de que 1 mm. Será 


di < 0,1 cm 


Suponhamos que o C.B.R. médio é de 10, 
e ter-se-á que 


2 
(Car) = Ss X01 


d(C.B.R.) «< 2 


Quer dizer, em pontos afastados de Im o 
C.B.R. não pode variar mais do que 2, isto é, 
de 20 “/o em relação à média. 

Mas se o C.B.R. for de 50 já as diferenças 


poderão ser de 


50º 
a >< 0,1 = S0 


quer dizer 100 º?'q em relação à média. A homo- 
geneidade terá portanto que ser tanto maior 
quanto pior for o terreno. 

Não tenho conhecimento de especificações que 
definam como esta homogeneidade há-de ser 
avaliada. No entanto a orientação que se vai 
definindo em Portugal (14) é a da aplicação ao 
problema dos métodos estatísticos de amostra- 
gem e ensaio. Os parâmetros a definir em cada 
caso serão: 


— Espaçamento da amostragem 

— Propriedades a determinar em cada amostra 

— Valores médius e máximas dispersões con- 
sentidas para tais propriedades. 


Enquanto porém tais especificações estatísti- 
cas não se encontrem definidas parece recomen- 
dável proceder sistemáticamente a uma homoge- 
neização dos materiais não só da base e sub-base 
como também da plataforma da estrada («sub- 
grade»). Esta homogeneização da plataforma é 
particularmente importante nas zonas de esca- 
vação, visto que nas zonas de aterro já houve 
alguma homogeneização nas operações de esca- 
vação, transporte, descarga e compactação. 
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A homogeneização pode consistir numa simples 
escavação, mistura e recompactação dos mate- 
riais no mesmo sítio. A profundidade a que deve 
ser levada esta escavação é muito variável de 
especificação para especificação. Nas diferentes 
regiões dos Estados Unidos da América vai-se 
desde 15 cm (6) a 75 cm (30) (4). 

A compactação exigida para esta plataforma 
é dum modo geral 95 “/o da compactação li- 
geira (AASHO T99); mas há no entanto esta- 
dos que exigem 100 “/, e mesmo 95 a 105 “/9, e 


outros que se contentam com 90 a 95 “/y (4). 


b) — Correcção granulométrica 


Para se conseguir que um material satisfaça a 
uma dada especificação granulométrica pode ser 
necessário corrigir a sua granulometria. Essa 
correcção poderá revestir-se das seguintes formas: 


1.º — Subtracção duma dada fracção granulo- 
métrica. 
2.º — Adição e mistura de outros materiais. 


A mera avaliação das granulometrias de mate- 
rial a subtrair ou a adicionar não oferece qual- 
quer dificuldade, uma vez conhecidas as curvas 
granulométricas dos diferentes materiais utilizá- 
veis. São problemas que se resolvem algebrica- 
mente ou graficamente com base nas curvas gra 
nulométricas ou na representação triangular. 

O que oferece dificuldades é projectar a ma- 
neira de efectuar a correcção de tal modo que o 
seu custo seja aceitável em face doutras possi- 


veis soluções do problema. Na sua generalidade 
o problema é muito complexo para poder ser 
abordado de maneira sistemática numa simples 
conferência. 

Resumidamente ele pode esquematizar-se con- 
forme se apresenta no Quadro VI. 


c) — Correcção da plasticidade 


Dados dois solos de índices de plasticidade Ki 
e Ke», respectivamente, em que proporção os 
deveremos misturar para obter um material com 
um dado índice de plasticidade K? 

Este problema pode ser aproximadamente re- 
solvido pelo método de Pimentel dos Santos (5) 
o qual consiste na aplicação da seguinte fórmula : 


— Ara Ki + As as ko 
Ara + As as 


K (10) 


na qual: 


— Ar e As = fracções em que os dois solos 
entram na mistura (Ai -H As == 1) 

— ai e az == constantes granulométricas dos 
solos a misturar definidas pela seguinte 
expressão: 

= 
à =— Ls (11) 
100 nm 
na qual 


— Y == ordenadas da curva granulométrica 
(percentagens do material que passa) cor- 
respondentes aos peneiros B. S. 200, 100, 52, 
25,14 e 7 ou os correspondentes americanos. 


QUADRO VI 


Correcções granulométricas de terrenos 


Pedras e calhaus ( => 60 mm) I 


Subtracção duma fracção gra- Seixos e areias 
nulométrica 


Areia fina, silte e argilas 


l 
| 


Seixos e areia 
Adição duma íracção granu- 
lométrica 


| Areia fina, silte e argila 
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| Crivagem mecânica ou escolha 
manual, seguida ou não de bri- 
tagem mecânica ou manual 


( Crivagem mecânica ou manual 


[| Por via seca mediante secagem, 

| desagregação e crivagem; ou 

| por via húmida mediante lava- 
gem e crivagem 


| Desagregação do terreno e mis- 
t tura do correctivo 


Desagregação do terreno e do 


t —Correctivo e respectiva mistura 


—n == número de ordenadas (no caso refe- 
rido = 6). 


A mesma fórmula permite também a previsão 
do limite de liquidez de mistura desde que Ky 
e Ks» representem as L. L. dos solos a misturar. 

Por este método podem portanto calcular-se as 
proporções em que se deverão misturar os solos 
disponíveis mas que tenham plasticidade em 
excesso (IP > 6) ou em defeito (IP 4) de 
modo a conseguir-se uma mistura que satisfaça 
as especificações. 

Modernamente utilizam-se outros processos de 
corrigir a plasticidade pela adição de cal ou 
cimento a que mais adiante me referirei. 


d) — Compactação 


Uma compactação pode ser definida especifi- 
cando o mínimo grau de compactação a obter 
depois de feita a compactação ou fixando a ma- 
neira como essa compactação deve ser feita 
(espessura das camadas, tipo e peso do cilindro, 
número de passagens, etc.). 


As especificações mais corréntes são do 1.º tipo; 
antigamente, porém, só havia especificações do 
2.º tipo. 

Qual deverá ser o grau de compactação a espe- 
cificar no projecto da obra? 

Como regra deverá especificar-se um grau de 
compactação que tire o maior partido da racional 
aplicação dos equipamentos usuais. 

Um dado equipamento, cujas principais carac- 
terísticas são o tipo de roda (roda lisa, pneus, 
pés de carneiro e o piso total), compactando 
camadas de espessura não superior a 20 cm rea- 
liza praticamente as suas possibilidades de com- 
pactação ao fim de cerca duma dezena de pas- 
sagens. Nessas condições, dado um terreno e um 
equipamento de compactação utilizado com efi- 
ciência e economia pode, em princípio, traçar-se 
no campo uma curva de compactação como a 
que se realiza no ensaio de laboratório. Haverá 
portanto no campo um peso específico aparente 
máximo e uma humidade óptima tal como no 
laboratório. A este resultado corresponde um 
certo custo. Querendo-se um grau de compac- 


QUADRO VII 


Comparação das máximas baridades secas (y) e dos teores óptimos de humidade (w) obtidos em ensaios 
de compactação de laboratório e de campo. Baseada em resultados obtidos pelo Road Research Laboratory 


Argila 


+ |wla 


(g/em?) | (9) | g/em?, 


a) — Ensaios de laboratório 


Compactação ligeira 

(AASHO normal) ....| 1,56 | 26. 
Compactação pesada 
(AASHO modificado)..| 1,80 | 17 | 
Relações de y e w (comp. | 
pes./comp. lig.) ...... 1,18 | 0,65 | 

b) — Ensaios de campo 

Cilindro de roda lisa 

(8 MOBO) paga pu secas 1.66 20 
Cilindro de pneus (12 

CORRE ssemas E ea dd RO 1,57 | 25 
Cilindro de pés de car- 

neiro (4,5 ton.) ...... 1a | 15 


Argila 
siltosa 


1,92 


1,18 


1,78 


1,66 


Areia : Seixo areno- 
; Areia h 
argilosa -argiloso 
w “ w Y w Y | w 


| (9/0) | (g/em?) (9/6) (g/em) 


21 1,84 14 | 1,94 11 | 2,06 9 

14 | 2,05 11 2,08 9 | 2,21 7 

0,67 | 1,12 | 0,78| 1,07 | 0,82| 1,07 | 0,78 

| 

16 1,86 14 2,11 8 2,21 7 

20 | 1,73 19 2,03 4) 2,02 7 

14 1,92 12 — — 2,05 5 
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tação bastante maior não é geralmente econó- 
mico diminuir a espessura das camadas nem au- 
mentar excessivamente o número de passagens. 
Terá portanto que se mudar de equipamento 
adoptando possivelmente outro mais pesado. 
O novo custo desta compactação será porém 
mais elevado. 

No quadro VII com base em ensaios feitos pelo 
Road Research Laboratory de Londres (6) com- 
param-se as características de compactação, má- 
xima baridade seca e teor óptimo de humidade, 
obtidas em cinco solos em ensaio de laboratório, 
de compactação ligeira e pesada, e em ensaio de 
campo com equipamento corrente, cilindro de 
roda lisa de 8 ton., cilindro de pneus de 10 ton. 
e cilindro de pés de carneiro de 4,5 ton. Por aí 
se vê que com tais equipamentos se obtêm graus 
de compactação à volta de 100 */o em relação à 
compactação ligeira e que portanto é razoável 
fazer essa exigência nos cadernos de encargos. 


QUADRO VII 


Comparação das máximas baridades secas (7) e dos 

teores óptimos de humidade (W) obtidos nos ensaios de 

laboratório e de campo. Baseada em resultados obtidos 
no Waterways Experiment Station 


Y | W 
(9) 


(g/cmº) 


a) — Ensaio de laboratório 


Compactação ligeira 
(AASHO normal) ...... 1,72 | 17,5 


Compactação intermédia 
(10 libras, 18”, 5 camadas, 
26 D./ camada) cu «us 1,79 | 16 


Compactação pesada 

(AASHO modificado) 

Relações de y e W com. 

pes./comp. lig......... 1,88 | 14 


b) — Ensaio de campo (16 passa- 
SORO spas o wa sa E 1,10 | 0,80 
Cilindro de pneus de 29 t 
5 EO) «ques a ssa 1,74 | 19 
Cilindro de pneus de 45 t 
(6,3: kelen) sau incas 178 | 17 
Cilindro de pneus de 57 t 
(10,5 KelonM usas o no: 1,88 | 15 
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Mas vê-se também que, com excepção dum caso 
(cilindro de roda lisa em seixo areno-argiloso), 
nunca se conseguiu obter 100 */o da compactação 
pesada. Aliás pelo mesmo quadro se vê que 
100 º/) da compactação pesada corrresponde a 
118 "/, da compactação normal para os terrenos 
argilosos, 112 "/y para a areia argilosa e 107 */o 
para as areias e seixos. 

Mudando de equipamento, por exemplo para 
cilindros de pneus de cerca de 29, 45 e 57 ton,, 
obteve a U. S. Corps of Engineers nos ensaios 
efectuados no Waterways Experiment Station (7) 
numa argila os resultados que se resumem no 
quadro VIII. Verifica-se aqui que somente se 
atingiu 100 */o da compactação pesada com o ci- 
lindro mais pesado. 

Estes dois exemplos dão uma ideia do que é 
razoável exigir dum dado equipamento e de 
quanto custa essa exigência. 


— €C. BR. 


Acabamos de fazer uma ideia de quanto custa 
subir as exigências de compactação. Para julgar 
da necesssidade dessa exigência torna-se neces- 
sário saber avaliar os benefícios que dela re- 
sultam. O método prático para atingir esse fim 
é por intermédio do C. B. R. Quanto maior é a 
compactação, maior será o €. B. R. e portanto 
menor a espessura do pavimento necessário. 
O custo desta diminuição de espessura do pavi- 
mento será portanto a contrapartida do excesso 
de custo da compactação. 

Para se avaliarem as vantagens das possíveis 
soluções torna-se necessário conhecer a influência 
no €. B. R. dum aumento do grau de compac- 
tação. Como essa relação varia de terreno para 
terreno torna-se necessário para cada caso fazer 
vários ensaios com diferentes energias de com- 
pactação, determinando o respectivo C. B. R. 
Na fig. 6 exemplifica-se a maneira de fazer essa 
avaliação com base em três energias de compac- 
tação segundo preconiza o U. S. Corps of Engi- 
neers (8). O €. B. R. pode avaliar-se por inter- 
polação. Assim para uma compactação de 95 º/% 
do AASHO modificado ter-se á um C. B. R. de 
cerca de 10º) enquanto que com uma compac- 
tação de 85 "/) AASHO modificado não se obterá 
mais do que um C. B.R. de 4º. A diferença 
de espessura de pavimento correspondente é de 


48 — 28== 20 cm, o que é como se vê muito 
importante. 
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12 lá 
TEOR DE HUMIDADE 


NOTAS : COMPACTAÇÃO: PILÃO 101. 
QUEDA 18º 
MOLDE CBR. 


EMBEBIÇÃO : 4 DIAS 


TEORES DE 
HUMIDADE DE 
COMPACTAÇÃO 


160 170 180 190 200 21 
BARIDADE SECA (9/em?) 


o —55 PANCADAS /CAMADA 

O-—25 PANCADAS /CAMADA 

A-— 10 PANCADAS /CAMADA 

O —TEOR DE HUMIDADE NO 
TOPO SUPERIOR DEPOIS 
DA EMBEBIÇÃO 


Fig. 6 — Exemplo da maneira de avaliar o C.B.R. do projecto 


3 — O caso das laterites 


Os materiais lateríticos ocorrem com grande 
frequência nas regiões tropicais e são por vezes 
os únicos materiais disponíveis em grandes ex- 
tensões. 

Nos materiais lateríticos distinguimos (2) as 
laterites em cascões ou concreções pisolíticas; os 
solos lateríticos, cuja fracção argilosa tem uma rela- 
ção molecular S5/0:/R:03<. 2 e os terrenos late- 
ríticos, nos quais há uma quantidade significativa 
de laterites ou de solos lateríticos. 

Na fig. 7 mostra-se a constituição química da 
fracção argilosa dos solos lateríticos. 

Estes solos lateríticos apresentam plasticidades 
relativamente elevadas e seriam regeitados em 


face da especificação americana. Verificou-se no 
entanto que mesmo com elevada plasticidade eles 
se comportam relativamente bem. Dada a impor- 
tância do problema para o Ultramar português, 
principalmente para Angola e Moçambique, fez-se 
um estudo do problema (2) visando o estabele- 
cimento de umas especificações. 

Verificou-se nesse estudo que tais solos apre- 
sentam uma pequena expansibilidade quando 
ensaiados no material que passa no peneiro N. 40, 
que é o mesmo material em que são feitos os 
limites de Atterberg. Na fig. 8, compara-se a 
sua expansibilidade com a de solos da região de 
Lisboa. Verifica-se que os solos lateríticos têm 
geralmente uma expansibilidade inferior a 5 “o, 
mesmo com 1. P. que para os solos de Lisboa cor- 
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O TERRENOS SEM CONCREÇÕES 


Fig. 7 — Constituição química da fracção argilosa dos solos lateríticos 


O TERRENOS LATERITCOS 


a 
AMA AS 


O TERRENOS NÃO LATERÍTICOS 
asa à TERRENOS NÃO LATERÍTICOS DE LISBOA 


EXPANSIBILIDADE “4 


Fig. 8 — Expansibilidade dos terrenos lateríticos comparada com a dos terrenos não lateríticos 


respondem a expansibilidades superiores a 10 º/o. Bases e sub-bases 

Modificaram-se por isso as respectivas especi- 
ficações e actualmente adopta-se, ainda a título L. E. < 40 (em vez de 25) 
experimental, a seguinte: LP S.15 (em vez de 6) 
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O NOVO TEODOLIO 


com calagem 
zenital automática 


O teodolito com os aperfei- 
çoamentos mais modernos, 
permitindo medições mais 
fáceis, rápidas e mais exactas 


| 


Automatismo duma concepção surpreen- 
dentemente simples: Prisma com líquido 
sem partes mecânicas, sem desgaste, ssm 
desarranjos, sem reparações. 


Detalhes no prospecto Th 154 E E: | 


WILD 


Modelo T1i-A C/imagem invertida 
Modelo T1i-AE C/imagem direita 


Solicitem catálogo ou demonstração 
ao representante exclusivo: 


WILD PORTUGAL, L.- 


LISBOA — PRAÇA DAS ÁGUAS LIVRES, 8, S/L 6 — Tel. 681127 


YBCNICA — XXI 


impresa Ultramarina de Sondagens & Fundaçãos, |. 


Associada: Empresa de Sondagens e Fundações 
TEIXEIRA DUARTE, LDA. 


SONDAGENS geológicas e geotécnicas — CAPTAÇÕES de águas 

subterrâneas — REBAIXAMENTOS do nível freático — CONSOLIDA- 

ÇÕES E ESTABILIZAÇÃO do solo — INJECÇÕES de cimento e 

outros produtos — INFRAESTRUTURAS de barragens, túneis e pontes 
eta FUNDAÇÕES de todos os tipos 


Única Empresa da especialidade com sede no Ultramar Português 


Sede-LUANDA 
Rua Serpa Pinto, 60 LISBOA L. MARQUES-C. P. 982 
| C.P. 2178 — Tel. 4226  P.da Figueira, 18, 3.º Esq. BEIRA — C. P. 1002 


| Companhia Portugusa de Trllara | -MPERMEÁRILIS“ 


S.A.R.L. 
ARMADURAS PRÉ-FABRICADAS MALHASOL O IMPERMEABILIZANTE OFICIALMENTE 
ACONSELHADO EM TODAS AS ARGA- 


| TETE MASSAS DE CIMENTO 
Ea Prédios, empenas, caves. — Depósitos para água, vinhos 
e apguardentes. — Facultam-se certificados de ensaio 
dos Laboratório Nacional de Engenharia Civil e da 

Junta Nacional do Vinho 


Tensão de segurança 3.000 Kg/cm? 


Distribuldor geral para todo o Império Português: 


AÇO Bl ANSELMO DE MATOS 
Av. Almirante Reis, 179, r/c. - Tel. 46439- LISBOA -1 


Revendedores em Lisbos: 


João Pereira Vareiro & Filhos, Lda. 
Rua Augusto José Vieira, 12 - Tel. 847480 


Manuel Pereira Matias, Lda. 
Tensão de segurança 4.000 Kg/cm? Rua de Pedrouços, r05-A - Tel. 6r 1153 


AÇOS DE ALTA RESISTÊNCIA PARA BETÃO ARMADO Manufacturas de Cimento Fortex, Lda. 


que pela eficiência técnica e económica a que conduzem, são 
um Indispensável elemento na moderna construção civil 
Para qualquer obra em betão armado consulte: 


DIAL — GABINETE TÉCNICO 


R. DO INSTITUTO INDUSTRIAL, 18-1.º DT.º — LISBOA 
TELEFONES 67 1224/5 E 673785 


Estrada de Benfica, 727 - Tel. 70 21 75 


Sequeira & Santos, Lda. 
R. José Joaquim Marques, 113 - Tel, 026105 - MONTIJO 
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Camadas de desgaste 


L. L. <, 40 (em vez de 35) 
I. P. < 6a 15 (emvez de 4a 9) 
Exp. € 5/p* 


Desgaste na máqnina de Los Ângeles 
< 65 “/ (em vez de 50º/) 


[HI— ESTABILIZAÇÃO COM CIMENTO 
1 — Especificações 


As especificações relativas ao solo-cimento ba- 
seiam-se, quanto ao solo, na granulometria, na plasti- 
cidade, na matéria orgânica e na existência doutras 
substâncias prejudiciais como os sulfatos, por exem- 
plo; quanto à dosagem em cimento baseiam-se 
as especificações na resistência à compressão e na 
durabilidade. 

Vejamo-las muito sumáriamente : 


a) — Granulometria 


O Road Research Laboratory (6) tendo em 
consideração a experiência americana indica os 
seguintes limites granulométricos para solos econôó- 
micamente estabilizáveis com cimento : 


— Dimensão máxima . .........css 75 mm (3º) 

— Passando no peneiro (B.S. 3/16' ou ASTM 4) 
4, 76 mm >> 50º 

— Passando no peneiro (B.S. 36 ou ASTM 40) 
0,42 mm >> 15º% 

— Passando no peneiro de (N.º 200) 0,074 mm 
< 509% 


Trata-se como se vê duma especificação bas- 
tante ampla cuja satisfação não apresenta geral- 
mente qualquer dificuldade. 


b) — Plasticidade 


O mesmo R.R.L. fixa os seguintes limites: 


L. L. << 40 
LP. <,:18 


que também são bastante amplos. No entanto o 
inconveniente dos solos mais plásticos será a 


* Aperfeiçoamentos recentes do ensaio de expansibi- 
lidade (Elda de Castro, L. N. E. C., Março de 1963), con- 
duziram à alteração deste valor para 10 */, (U. Nasci- 
mento, E. de Castro, M. Rodrigues, L. N. E. €., Março 
de 1963). 


maior quantidade de cimento necessária e a maior 
dificuldade na sua desagregação e pulverização 
para efeito de mistura do cimento. Peltier (9) 
refere os seguintes exemplos que são muito 
sugestivos: um solo com IP = 6 exigirá 8º/ de 
cimento, se IP==10 já necessitará de 10 º/% e se 
IP = 20 já carecerá de 15º/ o que, representando 
cerca de 300 kg/m”, já corresponde às dosagens 
do betão de cimento. 


c) — Matéria orgânica 


A influência da matéria orgânica é muito per- 
niciosa mesmo em pequena quantidade principal- 
mente por impedir a presa do cimento. São par- 
cularmente nocivos o açúcar e os ácidos húmicos, 
Em Moçambique, onde entre nós há maior expe- 
riência na utilização de solo-cimento, limita-se 
no respectivo Caderno de Encargos (10) o teor 
em matéria orgânica a 2º/. 

O tratamento prévio do terreno com cal, neu- 
tralizando os seus ácidos húmicos, pode tornar 
possível a sua ulterior utilização com cimento, 
Segundo indica Peltier (9) 5º/o de cal já é sufi- 
ciente para este efeito. 


d) — Sulfatos 


Os sulfatos são perigosos principalmente por- 
que atacam o cimento Portland normal. Em Mo- 
çambique (10) é interdita a utilização de solos 
cuja fracção argilosa contenha mais de 0,2º/o de 
sulfatos expressos em SO. 


e) — Resistência e durabilidade 


A resistência do“solo-cimento pode ser especi- 
ficada pela tensão de rotura à compressão simples 
ou pela Perda de peso em ensaios de durabilidade 


dade, originárias da Portland Cement Associa- 
tion (11). Na Inglaterra prevalecem as especifi- 
cações baseadas na resistência à compressão sim- 
ples. 

No quadro IX resumem-se as especificações 
referidas. 

As especificações americanas baseiam-se nos 
ensaios de compactação (ASTM — D 558-57 ou 
AASHO — T 134-57) de molhagem e secagem 
(ASTM — D 559-57 ou AASHO — T 135-57) e 
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TEN e PP 


QUADRO IX 
Especificações relativas ao solo-cimento 
Solos 
Ar, AY-4 | AZ.6, AZ 
Az-5e As | As, À; Ag, Às 
elo Siltosos | argilosos 
-arenosos 
Americanas (P. C. A.) | 
Durabilidade (7 dias) 
— Máxima perda de peso 
(12 ciclos) +. . «| 14% 10 “% 70% 
— Máxima variação de 
volume (12 ciclos) . - 2 º/ 2 º% 2º 


Resistência à compressão 


Solos arenosos 
14 a 21 kg/cm? 


Ensaio expedito (13) 
Mínima (7 dias) 


Conforme os ábacos respec- 
tivos (fig. 7) 
Resistência à compressão 
Inglesas (B. 5.) 


— Mínima (7 dias) . » » 
— Máxima (7 dias) +. . » 


17,6 kg/cm* 
25 kg cm* (9) 


Durabilidade 
— Resistência à compres- 
são após 7 dias de cura 
+. 7 dias de imersão, 
comparada com a re- 
sistência aos 14 dias 
de cura sem imersão. 
— Resistência após 14 ci- 
clos de gelo e desgelo 
comparada com a re- 
sistência de provetes 
padrão da mesma ida- 


dês mas Es 


> 80 9 (12) 


> 80 ºj9 (12) 


de gelo e desgelo (ASTM — D 560-57 ou AASHO 
— T 136-57). 

Como estes ensaios demoram mais de um 
mês a realizar a P.C. A. desenvolveu recente- 
mente um ensaio expedito (13) que tem por base 
a resistência à compressão simples dada por 
«ábacos». Estes ábacos foram traçados relacio- 
nando a resistência à compressão aos 7 dias de 
2229 solos com as especificações baseadas nos 
ensaios de durabilidade. 

O método só é válido para solos arenosos que 
contenham menos de 50 º/o de material inferior 
a 0,05 mm (silte mais argila) e menos de 20 “/y 
de material inferior a 0,05 mm (argila). 
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O método apresenta duas modalidades: o mé- 
todo A para solos que passam nó peneiro n.º 4 
(4,76 mm); e o método B para solos contendo 
material retido no peneiro n.º 4, 

A resistência mínima aos 7 dias tem que ser 
a indicada pelos ábacos que se indicam na fig. 7. 
Se for inferior terá que se efectuar o ensaio de 
durabilidade que fixará em definitivo a dosagem 
a adoptar. 

As especificações inglesas têm por base os 
ensaios da norma B.S. 1924:1957 (Methods of 
Test for Stabilized Soils). 


2 — Projecto duma estabilização com cimento 
a) — Dosagem com cimento 


A dosagem em cimento dum solo-cimento é 
feita ensaiando duas ou três dosagens, geral- 
mente afastadas de 2 “/o à roda da dosagem pro- 
vável para o tipo de solo em causa. Determina- 
das as resistências à compressão ou as perdas 
de peso no ensaio de durabilidade em função 
da dosagem, será fácil, por interpolação, calcu- 
lar a dosagem mínima correspondente à especi- 
ficação adoptada. 

Uma dosagem excessiva não convém para que 
o solo-cimento não funcione como laje formando 
portanto grandes fissuras de flexão e retracção. 
Peltier (9) indica 25 kg/cm? como limite supe- 
rior da resistência. 

As normas da P.C.A (13) indicam todos os 
passos a seguir para a determinação da dosa- 
gem mínima. 


b) — Dimensionamento 


O dimensionamento dos pavimentos em solo- 
-cimento costuma fazer-se pelo método C.B.R,, 
como os pavimentos flexíveis. No entanto, dada 
a maior rigidez do solo-cimento, há um certo 
efeito de laje, que melhora a distribuição das 
cargas no terreno de fundação. Há por isso uma 
certa tendência em tomar para espessura de cál- 
culo pelo método C. B. R. a espessura efectiva mul- 
tiplicada por um certo coeficiente. Estão em 
curso na Inglaterra experiências em troços expe- 
rimentais tendentes a fixar o valor deste coefi- 
ciente (12); entretanto, Peltier (9) sugere o valor 
médio de (4/3) quer dizer a espessura de cál- 
culo seria a espessura efectiva acrescida de 1/3. 


ência minima 
aos 7 dias, kg /cm? 


Resist 


Material inferior a 0,05 mm 
MÉTODO A — SOLOS PASSANDO NO PENEIRO N.º 4 (4,76 mm) 
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METODO B — SOLOS COM MATERIAL RETIDO NO PENEIRO N.º 4 (4,76 mm) 


Fig. 9 — Ábacos do método expedito da P, C. A. para fixação da resistência mínima 
aos 7 dias de solos arenosos estabilizados com cimento 


IV — ESTABILIZAÇÃO COM BETUME 


Os principais efeitos do betume numa estabi- 
lização de terrenos são o aumento da resistência, 
por efeito do aumento da coesão e a diminuição 
da susceptibilidade à água. 

A adição de betume pode ser feita por regas 
sucessivas de modo a obter-se uma penetração do 
material betuminoso na superfície do terreno, ou 
por mistura. No primeiro caso não se consegue 
estabilizar mais do que alguns milímetros de 
espessura, quando muito 1 a 2 centímetros, e 
pode consumir-se 1 a 5 |/m* de material betu- 
minoso; no segundo caso podem atingir-se es- 
pessuras de 10 a 20 cm com pleno controle da 
dosagem, assegurando-se um melhor aproveita- 
mento das características do material, 


É corrente dividir as estabilizações em dois ti- 
pos: o solo-betume, em que o terreno é um solo coe- 
sivo e a areia-betume em que o terreno é uma areia 
sem coesão. No primeiro tipo, o terreno já tem 
alguma coesão e, se se conservasse seco, não pre- 
cisaria provavelmente de qualquer estabilização ; 
a acção do betume é mais de defesa contra a água 
do que para aumentar a resistência mecânica. 

No segundo tipo, nas areias, cuja resistência 
mecânica é quase insensível à água, a acção do 
betume é dar-lhes coesão. 


1) — Especificações 


As especificações correntes dizem respeito ao 
solo, aos materiais betuminosos e à mistura. 
Vejamo-las sumáriamente : 
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a) — Granulometria e plasticidade 


Para os solo-betumes o Highway Research 
Board preconiza a seguinte especificação : 


4 (4,76 mm) > 50%% 
40 (0,42 mm) 35a100'/, 
» 200 (0,074mm) 10a50%/ 


Passando ao pen. n.º 


» » » » 


L. L.< 40 
É Fo 16 


b) — Material betuminoso 


Os materiais betuminosos utilizados são geral- 
mente os cut-backs e as emulsões. 

Os cut-backs, segundo as especificações ASTM, 
distinguem-se, como se sabe, pelas rapidez da 
cura, em «cut-backs» de cura rápida («rapid curing» 
ou R. C.) e de cura média («medium curing» ou 
M. C.) e de cura lenta («slow curing» ou S. C.) 
consoante o diluente betume é do tipo da gaso- 
lina, do petróleo ou do gasóleo, respectivamente. 
Cada um destes tipos tem várias graduações de 
viscosidade resultante quer da natureza do be- 
tume quer da sua concentração na diluição. 

A American Road Builders Association (A. R. 
B. A.) (15) deixa à escolha as graduações 1 a 4 
de qualquer dos tipos. No entanto para os solo- 
-betumes são mais correntes os MC. 1 e MC. 2, 

Para as emulsões a AASHO (M 140-49) espe- 
cifica três tipos de emulsões quanto à rotura: 
rápida (R. S.), média (M. S.) e lenta (S. S.) sub- 
dividindo os dois primeiros em várias gradua- 
ções de viscosidade e consistência. Para os mate- 
riais grosseiros, em grande parte retidos no 
peneiro de 3,18 mm (18”) e, práticamente sem 
material retido no peneiro 200 (0,074 mm), reco- 
menda as emulsões de rotura média (M. S.); 
para os materiais passando em grande parte no 
peneiro de 3,18 mm (1/8) e alguma quantidade 


no peneiro 200, recomenda a emulsão de rotura 
lenta S. 5. 


c) — Misturas 


As especificações para a mistura são geral- 
mente estabelecidas com base no ensaio Hubbard- 
-Fiel modificado e adoptado pela ASTM sob a 
designação D 915-47 T, 

Provetes com diâmetro e altura de 50,8 mm (2) 
são compactados sob 2722 kg (6000 libras' e 
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depois de submetidos a uma embebição em àgua 
durante 7 dias são submetidos a um ensaio de 
extrusão. 

As expecificações recomendam o seguinte: 


— Máxima absorção de água... 2º 

— Máximo aumento de volume 2º/% 

— Mínima resistência à extrusão: 
— com embebição 181 kg (400 libras) 
— sem embebição 454 kg (1000 libras) 


Há outras especificações que se baseiam no pene- 
trómetro de cone, preconizando uma resistêcia mi- 
nima de 80 kg/cm? (6), ou no ensaio de Cali- 
fornia que especifica uma resistência no mínimo 
de 1,75 kg cm? (25 ps). 


2) — Projecto duma mistura betuminosa 
a) — Dosagem 


Tal como o solo-cimento a escolha da dosa- 
gem a adoptar é feita ensaiando a absorção de 
água e a resistência de provetes preparados com 
várias dosagens. 


b) — Dimensionamento 


A espessura da camada é em geral calculada 
pelo método C. B. R. como pavimento flexível 
que é. 


V — OUTROS TIPOS DE ESTABILIZAÇÃO 
1 — Estabilização com cal 


A estabilização com cal, geralmente apagada e 
em pó, só por si ou misturada com pozolanas ou 
cinzas volantes está despertando muito interesse. 

A cal hidratada pode melhorar a plasticidade 
dos solos pouco plásticos, tais como o granito 
decomposto, seixo, etc., ou reduzir a plasticidade 
de argilas muito plásticas (16). Funciona portanto 
com um correctivo de plasticidade. 

Pode também ser utilizada como um aditivo para 
o solo-cimento fazendo aumentar a resistência 
deste e, por consequência, permitindo reduzir a 
sua dosagem (17). 

Pode ainda ser utilizada em substituição do 
cimento como agente estabilizador em percenta- 
gens da ordem de 6 “/o ou superiores. Tem no 
entanto o inconveniente de ser muito susceptível 
aos ciclos de gelo e desgelo (19 pág. 114). 
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Os mais modernos conhecimentos, obtidos na investigação 
científica e na experiência sob severas condições de serviço, 
também são considerados na construção dos nossos motores. 
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